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FOREWORD 


The  material  in  this  report  was  --limit ted  to  General  Research 
Corporation  by  Robert  Sampson,  a  GRC  c,  sultant,  and  is  based  upon  work 
which  he  and  T.  W.  Tuer  previously  perfc  .ned  at  the  University  of  Michigan 

Institute  of  Science  and  Technology. 

This  report  is  Volume  VI  of  the  Final  Report  for  the  Optical 
Signatures  Program  and  was  prepared  by  Thomas  M.  Zakrzewski,  Program 
Manager,  who  assumes  full  responsibility  for  the  revisions  of  the  contri¬ 
butions  herein.  The  program  is  under  the  technical  direction  of  W.  0. 
Davies  and  J.  L.  Hayes,  U.S.  Army  Advanced  Ballistic  Missile  Defense 
Agency,  Huntsville,  Alabama  office. 
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ABSTRACT 


The  exoatmospheric  target  signature  code  (EXOSIG) ,  a  segment  of 
the  Optical  Signatures  Code  (OSC)  ,  is  designed  to  calculate  the  optical 
signatures  of  targets  (reentry  vehicles,  decoys,  etc.)  during  exoatmos¬ 
pheric  flight.  EXOSIG  is  divided  into  nine  distinct  sections,  two  of 
which  supply  information  to  the  remaining  sections  and  not  to  the  user 
of  the  code.  These  two  routines  (TARG  and  EARTH)  are,  therefore,  pre¬ 
sented  as  appendixes  rather  than  in  the  main  body  of  this  volume.  The 
remainder  of  the  exoatmospheric  code  is  comprised  of  seven  primary  sec¬ 
tions,  each  of  which  is  available  to  the  user  depending  upon  which  types 
of  targets  and  calculations  he  desires.  Any  number  of  calls  may  be  made 
to  any  number  of  routines.  An  end-of-file  test  will  terminate  the  job. 
The  physics  utilized  in  the  development  of  each  section's  code  and 
the  input/output  information  necessary  to  or  resulting  from  each  of  the 
calculations  is  different  for  each  section  and  is  presented.  The  choice 
of  which  routine  (available  to  the  user)  is  required  by  him  is  based  upon 
the  target  configuration,  its  surface  material  properties,  and  whether 
the  user  desires  the  target's  optical  signature  or  the  polarization  con¬ 
tent  of  its  optical  signature. 
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I.  EXOATMOSPHERIC  TARGET  SIGNATURES  CODE  (EXOSIG)  USER'S  MANUAL 

The  exoatmospheric  target  signatures  code  (EXOSIG)  is  designed  to 
calculate  the  optical  signatures  of  targets  (reentry  vehicles,  decoys, 
etc.)  during  exoatmospheric  flight.  The  code  has  the  capability  of  cal¬ 
culating  the  heat  input  to  a  reentry  vehicle  during  launch  and  its 
resulting  optical  signature  during  launch  and  midcourse.  EXOSIG,  as  a 
whole,  handles  targets  reflecting  earthshine  specularly,  diffusely,  or 
bidirectionally  and  determines  both  the  self-emitted  and  reflected  con¬ 
tributions  to  the  optical  signature.  The  capability  for  calculating  the 
emitted  contribution  to  the  optical  signature  for  boosters  is  not  con¬ 
tained  in  this  code  because  of  the  uncertainty  in  the  heat  input  to  the 
booster  during  burn.  However,  the  reflected  contribution  to  the  optical 
signature  for  boosters  (often  times  the  most  significant)  can  be  calcu¬ 
lated  using  this  code. 

The  target  radiant  emission  at  a  specific  wavelength  or  in  a  wave¬ 
length  interval  is  determined  by  using  the  temperature  of  the  target  and 
assuming  a  Lambertian  radiator.  The  target  reflected  radiant  intensity 
at  a  specific  wavelength  or  in  a  particular  wavelength  interval  is  ob¬ 
tained  by  supplying  the  earth-emanating  radiation  at  the  particular 
wavelength  or  wavelength  interval  and  using  the  optical  material  proper¬ 
ties  for  the  particular  wavelength  or  wavelength  interval. 

The  exoatmospheric  code  has  the  capability  of  performing  seven 
separate  calculations,  therefore,  the  User's  Manual  is  divided  into 
seven  distinct  sections.  The  tasks  that  can  be  performed  by  using  part 
or  all  of  the  EXOSIG  routines  include: 

1.  The  calculation  of  the  temperature  and  the  resulting  radiant 
emission  from  a  reentry  vehicle  during  launch  amd  midcouise 
for  either  a  shrouded  or  unshrouded  launch  and  with  any 
launch  trajectory  (SURAD). 
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2.  The  calculations  of  the  contribution  to  the  optical  signature 
of  a  target  that  is  due  to  earth- emanating  radiation  that  is 
specularly  reflected  from  the  target  (SPEC1). 

3.  The  calculations  of  the  contribution  to  the  optical  signature 
of  a  target  that  is  due  to  earth-emanating  radiation  that  is 
diffusely  reflected  from  the  target  (DIFUS1), 

4.  The  calculations  of  the  equilibrium  temperature  for  a  target 
of  small  thermal  mass  and  the  resulting  directionally  emitted 
radiant  intensity.  In  addition,  the  contribution  to  the 
target  optical  signature  due  to  bidirectionally  reflected 
earth-emanated  radiation  is  calculated  (TINTF), 

5.  The  calculation  of  the  magnitude  of  the  error  introduced  in 
the  optical  signature  by  using  normal  material  properties 
instead  of  directional  material  properties  (DIREC) . 

6.  The  calculation  of  the  polarization  content  of  the  self-emitted 
radiant  intensity  (EP1). 

7.  The  calculation  of  the  polarization  content  of  the  radiant 
intensity  from  target  specularly  reflected  earth-emanating 
radiation  (SPOL) . 

Although  they  do  not  generate  significant  output  to  the  user,  two 
remaining  sections  (TARG  and  EARTH)  supply  information  essential  to  the 
operation  of  the  exoatmospheric  code  and  are  presented  as  Appendixes  I 
and  II.  The  relationship  of  TARG  and  EARTH  to  the  remainder  of  the  EXOSIG 
code  is  discussed  in  each  of  the  following  seven  sections.  In  addition, 
the  physics  utilized  in  the  development  of  each  section's  code  and  the 
INPUT/OUTPUT  information  necessary  to  or  resulting  from  each  of  the 
calculations  is  different  for  each  section  and  is  presented.  The  user 
must  decide  which  of  the  calculations  he  wishes  to  make.  This  decision 
will  be  based  upon  the  target  configuration,  its  surface  material  proper¬ 
ties,  and  whether  the  user  desires  the  target's  optical  signature  or 
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the  polarization  content  of  its  optical  signature.  The  user  should  now 
refer  to  the  section  of  this  manual  concerning  the  type  of  calculation 
he  desires. 
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11 •  REENTRY  VEHICLE  LAUNCH  AND  MIDCOURSE  TEMPERATURE  AND  RADIANT 

INTENSITY  CODE  (SURAD) 

A.  INTRODUCTION 

This  code  is  designed  to  calculate  the  surface  temperature  of 
a  reentry  vehicle  during  the  launch  and  midcourse  phases  of  ballistic 
flight.  This  surface  temperature  is  then  used  to  determine  the  radiant 
emission  from  a  reentry  vehicle  using  HRAD  (Appendix  III).  The  HRAD  code 
is  not,  at  present,  a  part  of  the  exoatmospher ic  code  and  must  be  called 
individually  through  the  OSC  driver  program.  The  general  procedure  for 
obtaining  the  surface  temperature  is  to  determine  the  flow  field  proper¬ 
ties  during  the  launch  phase  and  to  use  these  to  calculate  the  aerodynamic 
heat  input  to  the  reentry  vehicle.  The  net  heat  flux  is  obtained  by  a 
heat  balance  considering  aerodynamic,  solar  conduction  to  the  reentry 
vehicle  and  surface  radiation.  During  midcourse,  the  heat  balance  in¬ 
cludes  the  above  components  except  aerodynamic  heating  which  is  negligible. 
The  heat  shield  temperature  is  determined  using  the  net  heat  flux  by  a 
finite  difference  technique. 

The  procedure  allows  determination  of  the  surface  temperature 
for  any  body  shape  with  a  multilayered  heat  shield  on  any  arbitrary 
ballistic  trajectory.  Silo  blast  heating  and  heat  transfer  from  the 
booster  are  not  included;  however,  any  programmed  heat  input  can  be  con¬ 
sidered.  Ablation  and  nonequilibrium  effects  are  generally  not  appre¬ 
ciable  and  are  not  considered  in  this  analysis. 

B.  DEVELOPMENT  OF  EQUATIONS 

Any  description  of  how  launch  surface  temperature  is  calculated 
in  itself  is  a  description  of  the  entire  interaction  of  the  reentry 
vehicle  witu  its  environment.  One  must  begin  with  a  body  of  a  given  shape 
at  a  specific  velocity  and  altitude  since  this  defines  the  shock  condi¬ 
tions.  A  fluid  particle  of  air  must  be  followed  through  the  shock  and 
expanded  to  some  given  point  in  the  boundary  layer.  Heat  transfer  through 
fhe  boundary  layer  to  the  surface  must  then  be  calculated.  Midcourse 
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surface  temperature  is  much  easier  to  calculate  as  the  corrective  com¬ 
ponent  of  heat  transfer  is  zero  and  only  a  heat  balance  is  necessary. 

A  complete  flow  field  description  must  treat  regimes  from  free 
molecular  through  continuum  flow.  With  the  exception  of  a  narrow  transi¬ 
tional  regime  (see  Fig.  1),  the  flow  characteristics  are  calculable . ^ ^ 
Since  the  vehicle  spends  a  relatively  short  time  in  the  transitional 
regime,  transitional  theory  will  not  be  treated.  For  an  adequate  treat¬ 
ment  of  free  molecular  theory,  see  Ref.  1  or  2 . 

For  convenience,  the  discussion  of  continuum  flow  theory  is  pre¬ 
sented  for  blunt  vehicles  and  then  for  blunted  slender  vehicles.  For 
this  discussion,  the  blunt  body  is  defined  as  one  for  which  the  boundary 
layer  is  fed  by  near-normally  shocked  air  alone  (see  Fig.  2).  Thus,  the 
flow  field  of  interest  is  near  isentropic  and  surface  pressure  only  is 
needed  in  order  to  describe  the  expansion  from  the  stagnation  region  to 
any  given  point  in  the  boundary  layer.  Normal  shock  relations  for 
equilibrium  conditions  are  used  to  determine  properties  behind  the  shock. 

PcoVo,  “  P2V2  Continuity  (II-l) 


Momentum  (I 1-2) 


Energy  (II-3) 


P2  ’  »2RI2Z2 


State  (II-4) 


where  P  is  pressure 

V  is  velocity 

y  is  the  ratio  of  specific  heats  of  air 
M  is  Mach  number 
h  is  enthalpy 
R  is  the  gas  constant 
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T  is  temperature 
Z  is  the  compressibility 

2  (subscript)  denotes  properties  immediately  behind  the  shock 
00  (subscript)  denotes  ambient  properties 

3 

An  iterative  procedure  is  employed  to  determine  these  properties. 
Equilibrium  normal  shock  conditions  are  expanded  isentropically  to  a 
surface  pressure  determined  from  empirical  data  or  calculated  using  modi¬ 
fied  Newtonian  theory.  Although  modified  Newtonian  theory  breaks  down  at  the 
sphere— cone  juncture  for  small  cone  angles,  reasonable  accuracy  is  ob¬ 
tained  if  one  looks  at  distances  far  from  this  juncture.  The  cone-cone 
expansion  (if  the  vehicle  design  is  such)  is  similarly  treated.  State 
properties  for  equilibrium  air  are  those  published  by  Feldman.^ 

The  blunted- slender  body  is  one  for  which  the  boundary  layer  is 
is  fed  by  obliquely  shocked,  as  well  as  normally  shocked  air.  The  region 
of  transition  from  normally  shocked  air  to  obliquely  shocked  air  is  not 
considered,  only  points  before  and  after  transition.  Thus,  flow  prop¬ 
erties  at  the  external  side  of  the  boundary  layer  are  determined  before 

oblique  shock  influence  and  after  obliquely  shocked  air  dominates  boundary 
layer  flow. 

The  location  of  this  transition  region  is  determined  by  a  mass- 
flow-rate  balance  (see  Fig.  2). 


V  P«,VA  =  2"p  V  (<5  -  <5*)  r(x  ) 
el  1  11 


x2:  p«,VooA2  =  2ttp~  v~  -  6*)or(x„) 


e2  e2 


2  2 


(11-5) 
( II— 6) 


where  6  is  the  boundary-layer  thickness 

5*  is  the  displacement  thickness:  6*  =  f  !  1  -  — 2L_\  . 

0  \  peVe  J 
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pg  is  the  local  density  at  the  external  side  of  the  boundary 
layer 

Vg  is  the  local  velocity  at  the  external  side  of  the  boundary 
layer 

P^  is  the  free-stream  density 
is  the  free-stream  velocity 

r  is  the  radial  distance  from  the  longitudinal  axis  of  the 
reentry  vehicle 

Chapman"*  gives  6  and  6*  as  functions  of  surface  distance  x  for 

laminar  cones  at  a  variety  of  conditions.  and  are  the  areas 

defined  by  10°  deviations  from  the  normal  and  he  asymptotic  oblique 

shock,  respectively.  Shock  shapes  for  spheric  .1-nosed  bodies  are  obtained 

from  existing  experimental  and  theoretical  results.6  Thus,  the  limits  of 

this  transition  region  (x^  and  x^)  may  be  determined  from  the  above 

equations.  The  region  of  influence  of  normally  shocked  air  (x  <  x  ) 

7  ^ 

agrees  well  with  that  given  by  Chernyi. 

Since  we  know  the  limits  of  the  transition  region,  external 
flow  properties  for  the  regions  x  <  x^  and  x  >  may  be  determined  by 
expanding  isentropically  to  a  "Modified  Newtonian"  or  empirical  pressure 
from  normal  and  oblique  shock  properties,  respectively. 

The  calculation  of  the  surface  temperature  of  a  reentry  vehicle 
requires  determination  of  the  net  heat-transfer  rate  to  the  surface  and 
the  rate  of  conduction  of  the  heat  within  the  heat  shield.  The  net  heat- 
transfer  rate  to  the  surface  is  comprised  of  the  radiative  and  convective 
or  aerodynamic  contributions.  The  calculations  of  the  surface  radiation 
then  utilizes  the  calculated  surface  temperature  assuming  the  surface  to 
be  a  Lambertian  radiator  (Appendix  III). 

Semi-empirical  heat-transfer  relationships  are  available  for 
all  flow  regimes  with  the  exception  of  the  narrow  transitional  regime 
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discussed  in  inviscid  flow-field  theory.  However,  as  opposed  to  flow- 
field  theory,  heat-transfer  theory  may  be  used  in  its  full  form  to  perform 
inexpensive  calculations  on  simple  configurations.  Therefore,  approxi¬ 
mations  to  the  full  theory  will  not  be  discussed. 

Neglecting  the  thermal  motion  of  the  molecules,  the  convective 
heat-transfer  rate  in  the  free-molecular  flow  regime  can  be  expressed  as 


ccp  V  sin  0 


q  =  q 


conv 


where  a  is  the  accommodation  coefficient  (0  <  a  <  1) 
p  is  the  density 
V  is  the  velocity 
0  is  the  angle  of  attack 


(II-7) 


There  is  great  uncertainty  about  the  numerical  value  of  a.  No 

laboratory  tests  of  accommodation  coefficients  for  typical  reentry  vehicle 

materials  or  with  particle  energies  appropriate  to  reentry  velocities  have 

8  9 

been  published,  but  the  problem  is  discussed  by  Patterson,  and  Wachman. 
The  uncertainty  in  a  is  not  critical  because  the  body  is  in  appreciable 
free-molecular  heating  for  only  a  short  period  of  time. 


The  heat-transfer  rate  in  the  equilibrium  flow  region  of  the 
continuum  regime  is  needed  to  determine  the  heat  transfer  rate  for  the 
laminar  stagnation  point,  and  laminar  and  turbulent  flow  conditions  on 
other  portions  of  the  reentry  vehicle.  The  heat  transfer  to  the  stag¬ 
nation  point  as  given  by  Fay  and  Riddell,  is 


0.763,  ,0.1,  s0.4 

q-  =  ~o76(pij  ”  (pm)* 


Pr 


n„  D 

1  +  (Le°'52  -  1)— §- 


(h  -  h  ) 
s  w 


2(P  -  O 


1/2 
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where 


given 


where 


y  is  the  coefficient  of  viscosity  (lbf-sec-ft  ) 

Pr  is  the  Prandtl  number 
h  is  the  enthalpy  (BTU-lbm 
Le  is  the  Lewis  number 

h^  is  the  average  atomic-dissociation  energy  times  atom-mass 
fraction  in  external  flow 

is  the  noise  radius  (ft)  (radius  of  curvature  at  the  stagnation 
point) 

-2 

P  is  pressure  (lbf-ft  ) 
s  (subscript)  is  the  stagnation  point 
w  (subscript)  is  the  wall 

e  (subscript)  is  the  external  side  of  the  boundary  layer 


The  laminar  heat-transfer  rate  of  station  x  on  a  blunt  body,  as 
by  Avco-Everett,H  is 


*T  '  <Tw/T0)1/2(T0  +  202)/<Tw  +  2t>2) 

«  '  “o/RIo/X  W2 
0 

H  -  (2{/Ve)(dVe/d5)  -  |2«/ATPev2r2)^ve/dx) 

V2  * 

\  =  hs[4+he/hs(1-/L)]  =he+-f“ 
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x  =  wetted  length,  or  distance  transversed  by  following  the 
curved  path  of  the  body  surface  downstream  from  the 
stagnation  point  to  the  desired  station 

fi  =  recovery  factor 

r  =  radial  distance  from  the  longitudinal  axis  of  the  vehicle 


Although  derived  for  blunt  bodies,  (Eq.  II-9),  when  applied  to 

a  simple  cone,  gives  values  of  q  which  agree  with  values  obtained  from 

X  12 

an  empirical  relation  derived  by  AVCO-RAD. 


The  turbulent  heat-transfer  rate  at  station  x,  also  according  to 
AVCO-RAD,13  is 


-0.2.  .0.2 

x  (Rex) 


(1°g10Rex) 


2.58 


(11-10) 


in  which  Re  is  the  Reynolds  number  and  *  indicates  a  quantity  evaluated 

at  P  and  h*,  where 

e 

M.1  +  0.5g.1)  +  0.22,g.1) 

As  an  indication  of  the  similarity  of  heat-transf er  relationships  used  by 

investigators,  compare  the  following  relation  for  turbulent-heat-transfer 
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rate  as  used  by  GE,  with  that  of  AVCO-RAD  above: 


0.0296P  V  (y  ) 
e  e  e 


/  *  \0.  2/ 

lal)  [  £ 

if 

\V  /  \P 

\  e  /  ' 

j 

(r)0’  25 (h. 


h  ) 
w 


-’If 


„  0.25  1.25  . 

p  V  u  r  dx 
e  6  e 


k0.2 
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For  a  conical  body  (constant  flow  properties  on  cone),  this  reduces 


to 


,  (11-12) 

/  ft 

Further,  note  that  for  10  <  Re  <  10  ,  the  "Re  term"  in  Eq.  11-10  varies 

x  x 

from  0.156  to  0.189.  Thus,  Eq.  11-10  can  be  approximated  by 


(H-13) 

2 

where  (h  -  h  )  has  been  equated  to  V  /2.  The  AVCO  and  GE  equations  now 

S  0  6 

differ  by  a  constant  factor  of  only  1.14.  The  above  theory  adequately 
handles  both  blunt  and  slender  vehicles. 

Although  the  AVCO  equations  were  developed  for  hypersonic  flow, 

they  give  accurate  results  at  Mach  numbers  below  five.  At  Mach  three,  the 

AVCO  equations  agree  within  ten  percent  with  results  obtained  using  theorv 
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developed  by  E.  R.  Van  Driest  for  low  Mach  number  flight. 

The  previous  theoretical  treatments  have  assumed  equilibrium 
chemistry.  Nonequilibrium  effects  are  not  a  significant  factor  for  liquid 
fuel  booster  launches  and  are  significant  for  solid-fuel  booster  launches 
only  at  altitudes  after  peak  heating.  Therefore,  nonequilibrium  effects 
are  neglected  in  this  analysis. 

The  radiative  heat  transfer  consists  of  absorption  from  high- 
temperature  shocked  air  and  thermal  (surface)  emission.  Other  components, 
(solar  and  earth  emanating)  are  negligible  during  launch,  but  are  signi¬ 
ficant  during  midcourse. 
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In  the  stagnation  region,  the  maximum  radiative  heat-transfer 
from  high-temperature-shocked  air  encountered  by  current  reentry  vehicles 
is  approximately  one  order  of  magnitude  less  than  the  convective  heat- 
transfer.  At  positions  aft  of  the  stagnation  region,  the  ratio  of  radi¬ 
ative  to  convective  heat  transfer  is  much  lower  and  the  radiative  heat- 
transfer  may  be  neglected.  Therefore,  the  radiative  heat-transfer  is 
accounted  for  only  at  the  stagnation  point.  The  equation  used  is  an 
approximate  relationship  given  by  Detra  and  Hidalgo 


(11-14) 


2  -2 
where  K1  =  1.0  in  the  English  system  and  3.16  x  10  in  the  CGS  system; 

1  4  5 

and  K.  =  10  in  the  English  system  and  3.05  x  10  in  the  CGS  system.  A 
L  18 

similar  equation  is  given  by  Warzecha, 


q 


8 


OO 


where  C  is  a  constant  for  a  given  configuration  and  M  is  the  Mach 
number . 


Thermal  emission  is  extremely  important  in  maintaining  a  heat 
balance.  Indeed,  for  high-surface  temperature,  the  thermal  emission  q^ 
is  of  a  comparable  magnitude  to  the  aerodynamic  heating: 


.  A  C„/XT  \-1 

qe  -  C  e(X) A  5\e  W  If  dX 

0 

where  e(X)  is  the  spectral  emissivity  of  the  surface 
T^  is  the  wall  temperature 
X  is  the  wavelength 


(H-15) 
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4  -2  4 

C  =  3.7413  x  10  W-cm  -Mm 

C2  =  1.43  x  104  Mm  -  K 

If  we  combine  Eqs.  11-14  and  11-15,  the  net  heat-transfer  to  the  vehicle 
surface  during  launch  is  written 


CL/XT  \  1 

:  W  -  1/  dX  (11-16) 


and  q  is  given  by  Eqs.  II-7-10  depending  upon  the  flow  condition, 

aero 


During  midcourse  the  equation  for  heat  transfer  is 


net 


qearth  emanating  +  %olar 


e(X)X 


-5 


c9/xt 

2  w 


(11-17) 


qsolar  13  the  S°lar  radiati0n  conStant  or  zero*  However’  dearth  emanating 
varies  with  season,  time  of  day,  and  cloud  condition.  Since  the  reentry 

vehicle  has  a  large  mass  and  the  flight  time  is  of  relative  short  duration, 

the  effects  of  these  fluctuations  is  negligible  and  a  mean  value  is  assumed. 


Because  of  the  nonlinear  character  of  the  net  heat-transfer 
boundary  condition,  surface-temperature  determination  requires  a  computer 
solution  of  Fourier's  conduction  equation.  A  standard  finite  difference 
method  of  handling  this  problem  is  summarized  in  the  following  paragraph. 
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Since  a  steep  temperature  gradient  exists,  the  surface  temper¬ 
ature  must  be  calculated  by  considering  small  intervals  normal  to  the 
surface.  For  an  insulated  inner  surface,  the  temperature  at  all  intervals 
at  time  t  +  At  may  then  be  expressed  as  functions  of  the  heat  input  and 
the  interval  temperatures  at  time  t. 


The  first  interval  (outer  surface) 


(  q 

)  naero 

f  “rad 

1 

[■Wd  4y/2 

_ : 

^net  ^aero  qrad 
^net  ^cond  ^absorbed 


Transposing,  and  substituting  for  qabsorbed  and  qcon(J: 

l,t+At  A/2 


where 


A  =  C  p (Ay)2/kAt 
P 


K  is  the  thermal  conductivity  of  the  material 

C  is  the  specific  heat  of  the  material 
P 

p  is  the  density  of  the  material 


The  above  equation  expresses  the  unknown  surface  temperature  at  a 
future  time  (t  +  At)  as  a  function  of  the  variables 
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Tl,t+At  "  f^qnet*  Tl,t*  T2 ,t ^ 


The  heat  balance  must  be  formulated  in  subsequent  intervals  to  determine 
Tn,t+At* 

The  general  interval: 


n  -  1 


— 

1 

qcond 

qcond 

I 


T  ,+T  +  (A-  2)T 

n-1  n+1  n 


n,t+At 


(11-18) 


The  temperature  of  the  interval  of  a  material-bond  interface  is 
expressed  as  follows: 


n-1 


Ax 

2 

Ax 

2 


Material  1  _  n 

Bond  (Material  2) 


-n  +  1 


T 


n , t+At 


T  t  + 
n, 


T  -  T  \ 
n+l,t  n,t] 

Ax 


2At _ 

Ax  +  p„C  Ax 
1  Z  P2 


(11-19) 
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The  temperature  of  the  last  interval  of  the  structure  can  be  expressed  as 


T 


n , t+At 


(11-20) 


where  A  is  calculated  using  the  structure  properties. 

Temperatures  at  all  intervals  at  time  t+At  have  been  expressed 
as  functions  of  the  net  rate  of  heat  transfer  and  interval  temperatures 
at  time  t.  Since  the  net  heat  transfer  rate  is  a  function  of  surface 
temperature,  an  iterative  procedure  is  necessary  to  determine  surface 
temperature. 

C.  INPUT/OUTPUT  SPECIFICATIONS 

These  equations  are  combined  in  the  code  SURAD  for  calculating 
reentry  vehicle  launch  and  midcourse  surface  temperature.  Prior  to  using 
the  code,  the  reentry  vehicle  must  be  divided  into  isothermal  areas.  The 
temperature  for  a  central  point  (station)  in  each  area  will  be  calculated. 
The  heat  shield  at  each  station  is  divided  into  sublayers  with  the  thick¬ 
ness  of  each  sublayer  determined  by  the  requirement  that  A  (see  Eq.  11-17) 
be  greater  than  2.  This  is  necessary  when  using  a  finite  difference  tech¬ 
nique  for  temperature  calculation. 


Output  generated  by  SURAD  includes  a  time  versus  temperature 
array  of  each  layer  of  heat  shield  material  for  launch  and/or  midcourse. 
In  addition,  the  surface  temperature  from  the  code  can  be  put  into  HRAD 
(Appendix  III)  to  determine  the  radiant  emission  from  a  reentry  vehicle. 

The  input  data  necessary  for  using  this  code  is  tabulated  be¬ 
low  and  must  be  read  in  the  English  system  of  units.  As  many  data  sets 
as  the  user  requires  may  be  input  with  only  one  call  to  SURAD.  Since 
the  SURAD  code  runs  independently  of  all  remaining  routines  in  the  EXCSIG 
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package,  the  first  data  card  provided  for  input  to  SURAD  must  contain  the 
letters  SUR  in  the  first  three  columns  to  designate  that  the  operator 
desires  to  make  a  reentry  vehicle  surface  temperature  calculation. 

Since  arrays  LD  through  Z  defined  below  remain  constant  for  a 
majority  of  cases,  input  data  for  these  are  stored  on  a  permanent  file 
(See  Appendix  IV).  If  the  user  wishes  to  use  the  standard  input  to  arrays 
LD  through  Z,  he  must  set  the  "1TAPE"  input  variables  in  the  second  data 
card  to  "2".  If,  on  the  other  hand,  the  user  wishes  to  use  only  a  part 
or  none  of  the  standard  input,  he  must  set  the  "ITAPE"  input  variables  to 
"5"  for  each  of  the  arrays  LD  through  Z  which  is  to  use  nonstandard  input. 
A  definition  of  each  of  the  standard  arrays  is  given  below. 

Input  Data  F^r  SURAD 

LD(1) . . . LD(41)  are  forty-one  numbers  used  to  determine  the  density  of  the 
atmosphere  as  a  function  of  altitude  from  0  to  122  km. 

END3(1) . . .END3(44) ,  END4(1) . . .END4(44)  are  eighty-eight  numbers  which 

4 

represent  end  points  of  the  Mcllier  diagram  of  pressure  versus 
entropy.  These  end  points  are  necessary  to  indicate  when  to  use  the 
Mollier  diagram  for  determining  properties. 

LM(1) . . . LM(9 ) ,  ZPP (1) . . . ZPP (9) ,  TM(1) . . .TM(9)  are  three  fields,  each  of 

nine  numbers,  are  used  to  determine  the  temperature  of  the  atmosphere 
as  a  function  of  altitude  from  0  to  100  km. 

TB  is  the  lower  boundary  temperature  of  the  Mollier  diagram. 

R  is  the  universal  gas  constant. 

PO  is  reference  pressure. 

MINT42  and  EPSI42  are  the  increment  and  minimum  value  of  pressure  used 
in  the  tabulation  of  the  Mollier  diagram. 

Cl  and  C2  are  the  constants  of  the  first  degree  polynomial  for  enthalpy 
as  a  function  of  temperature. 
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ENDl ( 1 ) . . . END1 ( 36 ) ,  END2 (1) . . . END2 (36)  are  seventy-two  numbers  which  re¬ 
present  end  points  of  the  Mollier  diagram4  of  pressure  versus  entropy. 
These  end  points  are  necessary  to  indicate  when  to  use  the  Mollier 
diagram  for  determining  properties. 

RHOO  is  reference  density. 

MUO  is  reference  viscosity. 

RETRAN  is  the  theoretical  critical  Reynolds  Numbers. 

HB  is  the  base  enthalpy  of  the  Mollier  diagram. 

LE  is  Lewis  Number. 

PR  is  Prandtl  Number. 

TO  is  reference  temperature. 

IKL  is  the  number  of  iterations  allowed  on  the  hot  wall  heat  transfer 

loop  in  HTEMLP .  When  this  condition  is  exceeded,  the  iterative  loop 
is  completed  and  the  next  station  is  examined. 

GL  is  the  number  of  iterations  allowed  on  the  surface  temperature  loop  in 
HTEMLP.  When  this  condition  is  exceeded,  adequate  convergence  has 
not  been  obtained  and  the  program  is  terminated. 

ERROT  is  the  maximum  allowable  error  in  the  surface  temperature  iterative 
loop  in  HTEMLP. 

ERORT  is  the  maximum  allowable  error  in  the  hot  wall  heat  transfer  iterative 
loop  in  HTEMLP . 

TEMPER, H,Z  are  36  *  44  arrays  that  contain  the  temperature,  entalopy ,  and 
compressibility  factors  for  the  Mollier  diagram.  The  three 
tables  of  prints  are  discrete  points  on  a  Mollier  diagram  of 
pressure  versus  entropy. 

[The  remaining  data  will  vary  for  each  reentry  vehicle  and  launch  trajec¬ 
tory  and  must  be  supplied  by  the  user.] 
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IENO  is  an  end-of-data  test  parameter.  If  =  0,  data  follows.  If  =  9999, 
no  more  data  is  to  be  read  in. 

ERRORU  is  the  maximum  error  allowable  in  the  velocity  ratio  in  HALTLP 
in  percent. 

ACCOM  is  the  accommodation  coefficient. 

DT  is  the  time  increment  in  seconds  before  vehicle  left  atmosphere. 

FT  is  the  initial  flight  time  in  seconds. 

ENDFM  is  the  altitude  of  transition  from  continuum  to  free  molecular 
flow  in  kilofeet. 

ENDCON  is  the  altitude  of  the  end  of  launch  heating  and  the  beginning 
of  midcourse  flight  in  kilofeet.  * 

DE  is  an  array  of  cone  angles  for  oblique  shpck  in  degrees . 

KL  is  the  number  of  cone  angles  on  the  body. 

TIC  is  the  subscript  of  the  fir  altitude  that  launch  heating  is  to  be 
calculated  for. 

IL  Set  =*  2  If  both  sunlit  and  earth  heating  sections  of  midccurse  calcula¬ 
tions  are  desired. 

Set  =  1  If  only  sunlit  section  is  desired. 

Set  =  0  If  only  earth  heating  section  is  desired. 

LAUNCH  is  1  if  launch  heating  is  to  be  calculated  and  0  is  only  midcourse 
heating  is  to  be  calculated.  (1  =  unshrouded,  0  =  unshrouded) 

IOPRINT  is  the  printing  increment  for  all  cases  after  vehicle  has  left 
atmosphere. 

ICREMTIM  is  the  printing  increment  for  time  interval  before  vehicle  has 
left  atmosphere. 
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TL  is  the  subscript  of  the  last  altitude  and  velocity  to  be  input. 

TIR  is  the  subscript  of  the  first  velocity  altitude  input. 

KS(TIC) . . .KS(TL)  is  a  table  of  shock  angle  as  a  function  of  velocity  for 
a  specific  cone  angle  in  degrees. 

ALT (TIR) . . . ALT (TL) ,  VAM(TIR) . . . VAM(TL)  is  the  launch  altitude  and  velocity 
history  respectively  in  kilofeet  and  kilofeet/second . 

NL  is  the  number  of  stations  on  the  reentry  vehicle. 

NOBLI  is  the  front  station  feed  by  obliquely  shocked  air. 

RETURB  is  the  front  station  at  which  turbulent  flow  occurs. 

NSPHER  is  the  number  of  stations  on  the  sphere. 

NICONE  is  the  front  station  for  which  empirical  pressure  measurements  are 
available . 

RN  is  the  nose  radius  of  the  reentry  vehicles  in  feet. 

TRANSA  is  the  empirical  altitude  at  which  transition  to  laminar  flow  first 
occurs  in  kilofeet. 

TRNALT  is  the  empirical  altitude  at  which  transition  to  laminar  flow 
occurs  if  the  entire  reentry  vehicle  appears  to  become  laminar 
instantaneously  in  kilofeet. 

VAMOBL  is  the  velocity  at  which  the  flow  field  at  the  last  station  on  the 

reentry  vehicle  begins  to  be  fed  by  obliquely  shocked  air  in  kilofeet/ 
second . 

X  is  the  surface  distance  from  each  station  to  the  stagnation  point  in  feet 

XI  is  the  surface  distance  from  the  foremost  point  in  the  station  to  the 

stagnation  point  in  feet. 

DELTA  is  the  surface  angle  at  the  station  measured  from  the  axis  of 
symmetry  in  radians. 

RAD  is  the  radius  of  the  station  measured  from  the  axis  of  symmetry  in  feet 
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RI  is  the  radius  of  the  foremost  point  in  the  station  measured  from  the 
axis  of  symmetry  in  feet. 

TINI  is  the  initial  temperature  of  each  station  of  the  reentry  vehicle  in  R 
TSTL2  is  the  surface  temperature  at  (TIC-2)  in  °R. 

TITL2  is  the  second  layer  temperature  at  (TIC-2)  in  °R. 

FRAC  is  the  "gain"  in  the  wall  temperature.  Initially  it  is  equal  to 

1.0. 

2 

QH1  is  the  hot  wall  transfer  at  time  (TIC-1)  in  BTU/ft  -sec. 

NTI  is  the  time  of  the  first  calculation  in  midcourse  in  seconds. 

NTL  is  the  time  of  the  last  calculation  in  midcourse  in  seconds. 

PHI  is  the  absorptivity  of  the  outer  heat  shield  material  in  the  visible 
region. 

NDELT  is  the  time  increment  during  midcourse  in  seconds. 

DELSTR  is  the  thickness  of  the  structure  in  inches. 

3 

RHOSTR  is  the  density  of  the  structure  material  in  lbm/ft  . 

CPSTRU  is  the  heat  capacity  of  the  structure  material  in  BTlJ/lbm-°R. 

CONS'iR  is  che  thermal  conductivity  of  the  structure  material  in  BTU/ft  — 
°R-hr . 

DELY2  is  the  thickness  of  the  bond  between  the  first  and  second  layers 
of  heat  shield  material  in  inches. 

RH0M2  is  the  density  of  the  bond  material  between  the  first  and  second 

3 

layers  of  heat  shield  material  in  lbm/ft  . 

NTF  is  the  last  time  in  the  continuum  regime  before  entering  the  free 
molecular  regime.  It  is  the  time  of  the  last  altitude  before 
ENDFM  in  seconds. 

CP 2  is  the  heat  capacity  of  the  bond  material  between  the  first  and 
second  layers  of  heat  shield  material  in  BTU/lbm-  R. 
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C0N2  is  the  thermal  conductivity  of  the  bond  material  between  the  first 
and  second  layers  of  heat  shield  material  in  BTU/ft-  R~hr. 

DELY3  is  the  thickness  of  each  sublayer  of  the  second  layer  of  heat 
shield  material  in  inches. 

RH0M3  is  the  density  of  the  second  layer  of  heat  shield  material  in 
lbm/ft3. 

CP3  is  the  heat  capacity  of  the  second  layer  of  heat  shield  material  in 
BTU/lbm-°R. 

CON 3  is  the  thermal  conductivity  of  the  second  layer  of  heat  shield 
material  in  BTU/ft-°R-hr . 

DELY4  is  the  thickness  of  the  bond  between  the  second  and  third  layers 
of  heat  shield  materials  in  inches. 

RH0M4  is  the  density  of  the  bond  material  between  the  second  and  third 

3 

layers  of  heat  shield  material  in  lbm/ft  . 

CP 4  is  the  heat  capacity  of  the  bond  material  between  the  second  and 
third  layers  of  heat  shield  material  in  BTU/lbm-°R. 

C0N4  is  the  thermal  conductivity  of  the  bond  material  between  the  second 
and  third  layers  of  heat  shield  material  in  BTU/ft-  R-hr. 

DELY5  is  the  thickness  of  each  sublayer  of  the  third  layer  of  heat  shield 
material  in  inches. 

3 

RH0M5  is  the  density  of  the  third  layer  of  heat  shield  material  in  lbr /ft  . 

CP5  is  the  heat  capacity  of  the  third  layer  of  heat  shield  material  in 
BTU/lbm-°R. 

CON 5  is  the  thermal  conductivity  of  the  third  layer  of  heat  shield  material 
in  BTU/ft-0 R-hr. 

DELY6  is  the  thickness  of  the  bond  between  the  third  layer  of  heat  shield 
material  and  the  structure  in  inches. 

RH0M6  is  the  density  of  the  bond  between  the  third  layer  of  heat  shield 

3 

material  and  the  structure  in  lbm/ft  . 
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CP6  is  the  heat  capacity  of  the  bond  material  between  the  third  layer 
of  heat  shield  material  and  the  structure  in  bTU/lbm-°R.  j 

C0N6  is  the  thermal  conductivity  of  the  bond  material  between  the  third 
layer  of  heat  shield  material  and  the  structure  in  RTU/f t-°R-hr . 

ENDFL(l)  is  the  number  of  sublayers  of  the  fii.st  heat  shield  material. 

ENDTL  is  the  number  of  the  last  sublayer  of  the  second  heat  shield 
material. 

ENDFIF  is  the  number  of  the  last  sublayer  of  the  third  heat  shield 
material. 

NLL  is  the  total  number  of  layers  of  heat  shield  materia]. 

DELY  is  the  thickness  of  each  sublayer  of  the  outer  layer  of  heat  shield 
material  in  inches. 

RHOMAT  is  the  density  of  the  outer  layer  of  heat  shield  material  in 
lbm/ft3. 

CP  is  the  heat  capacity  of  the  outer  layer  of  heat  shield  material  in 
BTU/lbm-°R. 

CON  is  the  thermal  conductivity  of  the  outer  layer  of  heat  shield  material 
in  BTU/ft-°R-hr. 

EMIT  is  the  emissivity  of  the  outer  layer  of  heat  shield  material. 

TABLA  is  the  ablation  temperature  of  the  outer  layer  of  heat  shield 
material  in  °R. 

PRH  and  PRO  are  the  coefficients  of  the  linear  equation  used  to  calculate 
pressure  at  stations  NICONE  to  NL.  PRO  and  PRH  are  determined  from 
'  empirical  data  if  it  is  available. 
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111  *  CODE  FOR  EARTHSHINE  SPECULARLY  REFLECTED  FROM  TARGET  (SPEC1) 

A.  INTRODUCTION 

This  computer  code  is  designed  to  calculate  the  specularly  reflected 
earthshine  component  of  the  optical  radiation  signature  of  an  exoatmo- 
spheric  target.  The  general  procedure  is  to  subdivide  the  target  into 
elemental  areas  and  determine  the  contribution  from  each  element,  if  the 
earthshine  reflects  in  the  direction  of  the  viewer.  Input  tables  calcu¬ 
lated  by  TARG  (Appendix  I)  are  required,  so  the  assumptions  implicit  in 
that  code  are  also  pertinent  here.  In  addition,  the  target  is  assumed 
to  be  small  with  respect  to  its  distance  to  the  detector  and  to  the 
earth's  surface. 

B.  DEVELOPMENT  OF  EQUATIONS 

The  radiant  intensity  of  an  incremental  area  dA  of  the  target 
specularly  reflecting  radiance  L  (t}>)  from  an  earth  element,  toward  a 
detector  (Fig.  3)  can  be  written  as 

dl  =  pSL^ (4>) dA  cos  6  (III-l) 

g 

where  p  is  the  specular  reflectivity  of  the  material. 

To  find  the  total  radiant  intensity  due  to  reflected  earthshine, 
this  expression  is  integrated  over  all  of  the  target  that  is  visible 
to  the  detector  and  from  which  specular  reflection  from  the  earth  is 
possible.  The  angle  0  indicated  in  Fig.  3  is  determined  as  follows: 


R  =  V  - 


T 


(III-2) 


cos  9 


R  •  n/Rl 


(III-3) 


where  n  is  the  unit  normal  to  the  target  element,  R1  is  the  magnitude  of 

**►* 

R  and  is  obtained  from 
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3.  Specular  Reflection  Geometry 


UNCLASSIFIED 


UNCLASSIFIED 


R1  =  (R2  +  R2  +  R2)1/2  (III-4) 

x  y  z 

where  R  ,  R  ,  and  R  are  the  components  of  R  in  an  earth  fixed  coordinate 
x  y  z 

system. 

If  0  is  greater  than  90°,  there  is  no  radiation  to  the  viewer. 

A 

When  0  is  less  than  90°,  one  must  determine  the  vector  i  for  each  ele¬ 
ment  and  ensure  that  it  intersects  the  earth.  Three  independent  equations 

A 

are  required  in  order  to  determine  the  three  components  of  i.  Two 
independent  equations  are  determined  from  the  relationship  (Fig.  3) 

A  A  A  A 

n  x  r  =  i  x  n  (HI-5) 

and  the  third  from  the  relationship 

A  A  A  A 

r  •  n  =  n  •  i  (III-6) 

A 

To  ensure  that  the  vector  i  intersects  the  earth,  the  angle  0 
(Fig  3)  is  calculated  from  the  relationship 

->  * 

T  •  i  =  cos  0  (III-7) 

If  0  <  4),  where  4>  =  arc  sin  (earth  radius/earth  radius  and  target 

A 

altitude),  then  the  vector  i  intersects  the  earth  and  the  element  contri¬ 
butes  to  the  specularly  reflected  radiant  intensity  of  the  target. 

C.  INPUT /OUTPUT  SPECIFICATIONS 

Before  SPEC1  can  be  called,  TARG  must  be  requested  by  the  user 
(see  Appendix  I).  As  many  calls  as  required  may  then  be  made  to  SPEC1 
without  recalling  the  TARG  code. 

Table  1  defines  the  input  data  necessary  for  using  the  code.  The 
first  data  card  provided  for  input  to  SPEC1  must  contain  the  letters  SPE 
in  trie  first  three  columns  to  indicate  the  operator  desires  to  make  a 
specular  reflected  signature  calculation. 
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The  output  generated  by  SPEC1  includes  the  following: 

1.  Projected  area  (sq  mi) 

2.  View  factor  normalized  to  frontal  projected  area 

3.  View  factor  for  radiation  calculation  normalized  to  projected 
area 

4.  Nadir  angle  (degrees) 

5.  Aspect  angle  (degrees) 

6.  Radiant  intensity  (W/sr) 
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IV.  CODE  FOR  EARTHSHINE  DIFFUSELY  REFLECTED  FROM  TARGET  (DIFUS1) 

A.  INTRODUCTION 

This  computer  code  Ls  designed  to  calculate  the  diffusely  reflected 
earthshine  component  of  the  optical  radiation  signature  of  an  exoatmospheric 
target.  The  general  procedure  is  to  subdivide  the  target  and  the  portion 
of  the  earth  illuminating  it  into  elemental  areas.  A  double  numerical 
integration  is  then  performed  considering  the  elements  behave  as  flat 
plates.  Input  tables  calculated  by  other  codes  EARTH  (Appendix  II  and 
TARG  (Appendix  I)  are  required,  so  the  assumptions  implicit  in  them  are 
also  pertinent  here.  In  addition  the  target  is  assumed  to  be  small  with 
respect  to  its  distance  to  the  detector  and  to  the  earth's  surface. 

B.  DEVELOPMENT  OF  EQUATIONS 

The  radiant  intensity  of  an  incremental  area  dA  of  the  target 
diffusely  reflecting  incident  irradiance  dE  from  an  earth  element  and 
toward  a  detector  (see  Fig.  4)  can  be  written  as: 

dl  =  cos  0  cos  4>  (IV— 1) 

where  p  is  the  reflectance  of  the  material.  To  find  the  total  radiant 
intensity  (due  to  reflected  earthshine)  this  expression  is  integrated 
over  all  the  earth  which  is  illuminating  the  target  and  over  all  the  target 
which  is  visible  to  the  detector.  Using  vector  algebra  the  cosine  terms 
are  simply 


cos  9  =  n  •  i  (IV-2) 

cos  <f>  =  n  r  (IV-3) 
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Figure  4.  Reflection  Geometry 
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~r  ,*r  i  * 

where  n  is  the  unit  normal  vector  of  the  element,  (n  =  A/|A|),  i  is  a 
unit  vector  toward  the  source  of  the  incident  radiation  (e.g.,  an  earth 
element),  and  r  is  a  unit  vector  toward  the  detector. 

Given  the  latitude  ip,  longitude  X,  and  altitude  6  (normalized  to 
earth  radius)  of  the  target  (subscript  t)  and  the  detector,  (subscript  d) 
their  position  vectors  and  may  be  written  in  an  earth-fixed 
coordinate  system  as: 


•f 

R  -  0  fe  cos 
t  t  x 

-y 

R  =  Me  cos 
d  d  x 


cos 

Ip,  cos 
a 


X  +  e  cos 
t  y  t 


X  ,  +  e  cos  \p. 
ay  d 


sin  X  +  e  sin  ip.  )  (IV-4) 
t  z  t 

sin  X  +  e  sin  tK)  (IV- 5) 
d  z  a 


where  e  ,  e  ,  and  e  are  shown  in  Fig.  5.  Using  these  position  vectors, 
x  y  z 

the  vector  r  is  simply: 

r  =  R  -  R .  (IV-6) 

t  a 

and 

r=r/|r|  (IV-7) 

Alternatively,  if  the  nadir  angle  5  and  the  azimuth  angle  x  to  the 
viewer  relative  to  the  target  is  specified  (see  Fig.  6)  we  can  write 
simply: 


r  = 


e  sin  6 

x 


cos 


X  +  e^  sin  6  sin  x 


§  cos  6 
z 


(IV-8) 


where  the  target  is  located  on  the  z  axis. 
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A 


Figure  5.  Earth-Fixed  Coordinate  System 
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'1 


Figure  6.  Alternative  Earth-Fixed 
Coordinate  System 
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As  a  simple  check  on  the  code,  the  projected  area  A  of  the  target 

P 

is  also  calculated.  This  may  be  written  as: 


cos  <f>  dA 


(IV-9) 


For  computer  evaluation,  qs.  IV-1  and  IV-9  are  converted  to  finite 
difference  form: 


1  =  *  S  AAj  COS  *j  2  AEk  C0S  9jk  (IV-10) 

J  k 

and 

A  =  ^  AA  cos  <t>  ■  (IV-11) 

J 

where  the  only  elements  considered  are  those  for  which  the  cosine  terms 
are  positive. 

C.  INPUT/OUTPUT  SPECIFICATIONS 

Before  DIFUS1  can  be  called,  TARG  and  EARTH  (in  that  order)  must  be 
requested  by  the  user  (see  Appendixes  I  and  II).  As  many  calls  as  required 
may  then  be  made  to  DIFUS1  without  recalling  the  TARG  and  EARTH  codes. 

Table  2  defines  the  input  data  necessary  for  using  this  code.  The  first 
data  card  provided  for  input  to  DIFUS1  must  contain  the  letters  DIF  in  the 
first  three  columns  to  indicate  the  operator  desires  to  make  a  diffusely 
reflected  signature  calculation. 

The  output  generated  by  DIFUS1  includes  the  following: 

1.  Projected  area  (sq  mi) 

2.  View  factor  normalized  to  frontal  projected  area 

3.  View  factor  for  radiation  calculation  normalized  to  projected  area 

4.  Nadir  angle  (degrees) 

5.  Aspect  angle  (degrees) 

6.  Azimuth  angle  (degrees) 

7.  Radiant  intensity  (W/sr) 
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V.  CODE  FOR  DIRECTIONALLY  EMITTED  AND  BIDIRECTIONALLY  REFLECTED 

RADIATION  FROM  TARGETS  (TINTF) 

A.  INTRODUCTION 

This  computer  code  is  designed  to  calculate  the  directionally 
emitted  raiiation  plus  the  bidirectionally  reflected  earthshine  component 
of  the  optical  radiation  signature  of  an  exoatmospheric  target.  The 
general  procedure  is  to  subdivide  the  target  and  the  portion  of  the  earth 
illuminating  it  into  elemental  areas.  The  emitted  radiant  intensity  in 
the  direction  of  the  viewer  is  calculated  using  input  directional  emissi- 
vity.  The  reflected  earthshine  in  the  direction  of  the  viewer  is  cal¬ 
culated  using  a  double  numerical  integration  considering  the  earth  and 
target  elements  behave  as  flat  plates.  The  bidirectional  reflectivity 
is  calculated  using  a  semi-empirical  technique  that  requires  the  refrac¬ 
tive  index  absorption  coefficient,  and  normal  diffuse  reflectance  for 
the  material  as  well  as  a  measured  distribution  function  for  specular 
reflectance.  Input  tables  calculated  by  other  codes  EARTH  (Appendix  II) 
and  TARC  (Appendix  I)  are  required,  so  the  assumptions  implicit  in  them 
are  also  pertinent  here. 

B.  DEVELOPMENT  OF  EQUATIONS 

The  radiant  intensity  of  an  incremental  area  of  dA  of  the  target 
emitted  toward  a  detector  can  be  written  as: 

dl  =  e(0)  dA  cos  0  (V-l) 

where  e  (0)  is  the  directional  emittance  of  the  material.  To  find  the 
total  emitted  radiant  intensity,  this  expression  is  integrated  over  all 
the  target  which  is  visible  to  the  detector.  The  angle  0  indicated  in 
Fig.  3  is  determined  as  follows: 

3  =  V  -  ^  (V-2) 

0  =  t  •  n/Rl 


42 

UNCLASSIFIED 


cos 


(V-3) 


UNCLASSIFIED 


where  n  is  the  normal  to  the  target  element  R1  is  the  magnitude  of  ^ 
and  is  obtained  from 

R1  =  \/r2  +  K2  +  R2  (V-4) 

V  x  y  z 

where  R  R  and  R  are  the  components  of  R  in  an  earth  fixed  coordinate 
x  y  z 

system. 

If  6  is  greater  than  90°,  there  r'.s  no  radiation  to  the  viewer.  The 
program  also  checks  for  earth  blockage  of  the  radiation  from  the  target  to 
the  viewer  by  calculating  the  angle  between  a  vector  from  the  viewer  to 

the  center  of  the  earth  and  from  the  target  element  to  the  viewer.  This 

-y 

angle  is  used  to  calculate  the  length  of  a  normal  vector  to  the  R  vector 
to  the  center  of  the  earth.  If  either  4>  is  greater  than  90  or  K  is 
greater  than  the  radius  of  the  earth,  or  if  6  is  less  than  u,  there  is  no 
earth  blockage  (Fig.  7). 

The  radiant  intensity  of  an  incremental  area  dA  of  the  target  bi¬ 
directionally  reflecting  irradiance  dE  from  an  earth  element,  toward  a 
detector  (Fig.  8)  can  be  written  as: 

dl  =  p(<P,  9,  <t>)  dE  cos  tjj  dA  cos  6  (V-5) 

where  p(<J>,  0,  <f>)  is  the  directional  reflectivity.  To  find  the  total 
reflected  radiant  intensity,  this  expression  is  integrated  over  all  the 
earth  that  is  illuminating  the  target  and  over  all  the  target  that  is 
visible  to  the  detector.  The  angle  9  has  been  determined  previously. 

The  angle  ip  is  obtained  from: 

n  •  e  =  cos  ip 

The  angle  <p  is  obtained  as  follows: 
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Figure  8.  Bidirectional  Reflectance 
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A  vector  d  is  determined  by  the  relationship  (see  Fig.  8)  : 

r  x  n  =  d  (V-7) 

A  vector  c  is  determined  from  the  relationship  (see  Fig.  8): 


e  x  n  =  c 


(V-8) 


The  angle  <p  is  then  obtained  from  the  relationship! 


cos  <f>  =  c  •  d 


(V-9) 


The  reflectivity  is  obtained  from  the  following  relationship: 


p(e,  \p,  <t>)  = 


s  2 

P  (0  )  *  cos  8  * 

n 


_  8  ,  S 
PL  +  PP 


n  2  9 

2  cos  i p  cos  0  (n  -  1)  +  k 


(n  +  l)2  +  <2 

(V-10) 


g 

where  p  (0  ) 


n 

K 

d 

P 


is  the  directional  specular  reflectance  (0  is  indicated 

n 

in  Fig.  8) 

are  the  Fresnel  reflectances  calculated  using  8  in  Fig.  8 
is  the  material  refractive  index 
is  the  material  absorption  coefficient 
is  the  diffuse  reflectance  of  the  material 


C.  INPUT/OUTPUT  SPECIFICATION 

Before  TINTF  can  be  called,  TARG  and  EARTH  (in  that  order)  must  be 
requested  by  the  user  (see  Appendixes  I  and  II).  As  many  calls  as  re¬ 
quired  may  then  be  made  to  TINTF  without  recalling  the  TARG  and  EARTH 
codes. 
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Table  3  defines  the  necessary  input  data  for  using  this  code.  The 
first  data  card  provided  for  input  to  TINTF  roust  contain  the  letters  TIN 
in  the  first  three  columns  to  indicate  the  operator  desires  to  make  a 
calculation  using  this  code.  The  output  generated  by  TINTF  follows: 

1.  Projected  area  of  target  as  seen  by  viewer  (sq  mi) 

2.  Area  of  target  as  seen  by  both  earth  and  viewer  (sq  mi) 

3.  Radiant  intensity  of  target  toward  the  viewer  due  to  reflected 
earthshine,  assuming  a  bidirectional  reflectance  of  1  (W/sr) 

4.  Quantity  needed  to  determine  target  emitted  radiant  intensity 
to  the  viewer 

5.  Bidirectional  reflected  earthshine  radiant  intensity  in  the 
direction  of  the  viewer  (W/sr) 

6.  Bidirectional  target-reflected  earth-reflected  sunlit  radiant 
intensity  in  direction  of  the  viewer  (W/sr) 

7.  Total  overall  wavelengths  of  emitted  radiant  intensity  of  objects 
in  direction  of  viewer  (W/sr) 

8.  Reflected  radiant  intensity  in  direction  of  viewer  (W/sr) 

9.  Irradiance  at  target  due  to  earthshine 

10.  Irradiance  at  target  due  to  reflected  sunlight 

2 

11.  Sigma  temperature  to  the  fourth  power  (W/sr-cm  ) 

12.  Sum  of  total  emitted  and  reflected  radiant  intensity  (this  is 
useful  only  when  reflected  signature  is  a  total  overall  wave¬ 
length) 

13.  Temperature  (K) 

If  only  target-emitted  radiant  intensity  is  desired,  set  NOREFL 
equal  to  1,  and  only  target  equilibirum  temperature  and  effective  emissi- 
vity  in  the  direction  of  the  viewer  are  calculated.  To  obtain  the  emitted 
radiant  intensity  in  a  specific  wavelength  interval,  the  output  tempera¬ 
ture  and  effective  emissivity  must  be  put  into  HRAD  (Appendix  III)  and 
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the  result  then  added  to  the  reflected  component.  The  HRAD  code  is  not 
at  present,  a  part  of  the  exoatmospheric  code  and  must  be  called  indivi 
dually  through  the  OSC  driver  program. 


48 

UNCLASSIFIED 


TABLE  3 

INPUT  DATA  FOR  TINTF 
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TABLE  3  (Cont.) 
INPUT  DATA  FOR  TINTF 
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TABLE  3  (Cont.) 
INPUT  DATA  FOR  TINTF 


Values  Used 
if 

Not  Set 

0.14 

o 

H 

tu 

H-J  tH 

0^1  01 
(0  tu 
a)  >  0 

< 

• 

c 

60  O  r-H 

13  > 

• 

• 

S3 

CM 

CN 

B 

e 

CO 

u 

n 

n 

•H 

CO 

CO 

d 

1 — ^ 

tu 

3 

3 

c 

o 

C 

a 

•H 

o 

0 

U 

•H 

•H 

•H 

•u 

■U 

c 

Cd 

cd 

•H 

•H 

M-l 

T3 

TJ 

2 

cd 

U 

Cd 

u 

JC 

cd 

U 

u 

O 

cd 

in 

w 

SIFIED 


UNCLASSIFIED 


VI.  DIRECTIONAL  PROPERTIES  EFFECT  CODE  (DIREC) 

A.  INTRODUCTION 

Computer  code  DIREC  was  designed  to  determine  the  effect  of 
directionality  in  the  optical  properties  on  the  optical  signatures  of 
certain  targets.  These  signatures  are  comprised  of  both  self  emitted 
and  reflected  earthshine  so  that  both  the  emissivity  and  the  reflectivity 
are  required.  Fresnel  equations  are  employed  to  evaluate  the  optical 
properties.  The  calculated  parameters  from  this  code  are  total  effective 
emissivity  and  reflectivity  as  a  function  of  the  viewer-target-earth  geom¬ 
etry.  Two  simplifying  assumptions  were  made  in  this  code:  uniform 
temperature  and  optical  properties  over  the  surface  of  the  target  and  uni¬ 
form  earth  radiance  distribution.  Either  of  these  could  be  relaxed  with 
little  conceptual  difficulty  but  the  complexity  of  the  code  would  be 
increased  considerably. 

Options  are  provided  for  inputting  either  the  viewer  position  (in 
latitude,  longitude,  and  altitude),  or  the  viewer  nadir  and  azimuth  angles 
relative  to  the  target.  Results  of  another  code  TARG  (Appendix  I)  are 
required,  so  the  assumptions  implicit  in  it  are  also  pertinent  here. 

B.  DEVELOPMENT  OF  EQUATIONS 

If  the  viewer  latitude  0|iy)  ,  longitude  (Ay)  and  normalized  altitude 
(By)  are  specified,  the  components  of  the  viewer  position  vector  V  are 
found  fvutn: 


x  =  8tj  cos  cos  A 
V  V  V  v 

yy  =  Sy  cos  ^y  sin  Ay 


z 


V 


gy  sin  i|iy 
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These  components  are  referenced  to  an  earth-fixed,  right-handed  coordi¬ 
nate  system  with  the  origin  at  center  of  the  earth,  the  z-axis  toward 
the  North  Pole  and  the  x-axis  toward  the  prime  meridian  (see  Fig.  9). 
Similarly  for  the  target  position  vector  T: 

XT  =  6T  COS  COS  XT 
yT  =  6t  cos  4»t  sin  Ax 

zT  =  6t  sin  i|^T 

The  range  vector  R  from  the  target  to  the  viewer  is  simply  the  difference 
of  these  two: 

R  =  V  -  T 

For  the  option  of  specifying  the  viewer  nadir  6  and  azimuth  X 
angles,  the  target  must  be  located  along  the  z-axis  and  the  target's  axis 
in  the  x-z  plane  (see  Fig.  10).  Here  different  equations  are  used  to 
determine  R: 


R  =  sin  6  cos  X 
x 


R  =  sin  6  sin  X 

y 


R  =  -cos  6 
z 


TARG  (Appendix  I)  provides  area  normal  vector  (A)  components  from 
each  target  element  in  the  same  coordinate  frame.  The  area  an  element 
projects  toward  the  viewer  (A  )  is  (see  Fig.  11): 
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A 


Figure  9.  Viewer- 
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Figure  11. 


Reflection  Geometry  for  an  Elemental  Area 
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A  =  A  •  R 
P 

where  R  denotes  a  unitized  vector  obtained  from: 

R  =  R/  |  R  | 

If  this  projected  area  is  negative,  the  element  faces  away  from  the 
viewer  and  thus  is  neglected.  The  angle  of  reflection  (0  ) ,  may  be  found 
from: 


cos  0  =  A  /  A 

r  P 


19 


The  Fresnel  equations  give  the  reflectivities  ^  ’  and  Pp 
inducing  parallel  "L"  and  perpendicular  "P"  polarization: 


2  2  2 

(n  cos  0^  -  cos  x)  +  <  cos  0^ 

Pf  =  '  2  2 

(n  cos  0^  +  cos  x)  +  K  cos  0r 


\  2  2  2 

(n  cos  x  “  cos  0r)  +  K  cos  X 

Pp  =  ~  2  2 

(n  cos  x  +  cos  9  )  +  K  cos  X 


where  n  is  the  refractive  index  and  k  the  extinction  coefficient  of  the 
material,  and: 


sin  x  =  “  sin  0 
A  n  r 


The  unpolarized  reflectivity  p  is  simply  the  average  of  those: 


P  = 
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and  the  emissivity  "e"  (for  an  opaque  material): 
e  =  1  -  p 


Assuming  the  target's  surface  to  be  at  a  uniform  temperature,  the 
effective  emissivity  e  is  simply: 


e  = 


N 

£ 

i=l 


e  A 
1  P 


P(A  ) 


where  N  is  the  total  number  of  elements,  P  is  a  function  which  is  unity 

if  its  argument  is  positive  and  zero  if  nonpositive,  and  the  total  pro¬ 
jected  area  A  is: 


*  £  A„  p«  ) 


i=  1 


The  effective  reflectivity  p  is  similarly  expressed  assuming  the  earth- 
shine  to  have  a  uniform  spatial  distribution: 


P  = 


l  N 

r£ 

r  i=l 


P.A  P(A  )P(3!) 
1  Pi  Pi  1 


where  the  total  reflecting  area  A 


is  defined  by: 
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The  function  P^)  indicates  whether  the  ith  element  actually  reflects 
the  earth's  image  to  the  viewer.  The  argument  3'  is  defined: 

3'  -  <p  -  3 

m 

where  ^  is  the  nadir  angle  to  the  earth's  horizon  and  3  is  the  angle 
between  the  nadir  position  and  a  unit  vector  I  indicating  the  incidence 
direction  (see  Fig.  12). 

They  are  evaluated  from 

cos  3  =  -I  •  T 


and 


R 

e 


where  is  the  radius  of  the  earth  and  Ht  is  the  altitude  of  the  target 
from  the  surface  of  the  earth. 

Since  for  total  reflection,  the  angle  of  incidence  0^^  equals  the 

angle  of  reflection  9  : 

r 

A  x  I  =  R  x  A 


and 


I  =  cos  8, 
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PORTION  OF  EARTH 
VISIBLE  FROM  TARGE 


TARGET 


Figure  12.  Reflection  Geometry  (Nadir  Angle) 
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After  expanding  we  get  four  equations  (three  are  dependent)  in  the  three 
unknown  components  of  I,  which  may  be  solved  to  give: 


I  =  A  cos  0  +  C.A 
XX  2  z 

Ty  =  Ay  COS  0  "  CxAz 

I  =  A  cos  0  +  C,A 
z  z  1  y 


+  C.A 
3  y 


-  C.A 
3  x 


-  C  A 
2  x 


where  C  =  A  R  -  A  R 
1  z  y  y  z 


C  =  A  R  -  A  R 
2  x  z  z  x 


C3  5  <C1A*  +  Vy>/A2 


or  if  A  =  0, 
z 


=  A  R  -  A  R 
3  z  y  y  x 

Finally,  the  aspect  angle  0  (measured  from  nose  on)  is  found  from  the 
vector  specifying  the  target's  axis  N  and  the  range  vector  R: 

cos  0  =  N  •  R 

C.  INPUT/OUTPUT  SPECIFICATIONS 

Before  DIREC  can  be  called,  TARG  must  be  requested  by  the  user  (see 
Appendix  I).  As  many  calls  as  required  may  then  be  made  to  DIREC  without 
recalling  the  TARG  code.  Table  4  defines  the  necessary  input  for  using 
this  code.  The  first  data  card  provided  for  input  to  DIREC  must  contain 
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the  letters  DIR  in  the  first  three  columns  to  indicate  the  operator  desires 
to  use  this  code.  The  output  generated  by  DIREC  follows: 

1.  Projected  area  (sq  mi) 

2.  Reflected  area  to  projected  area  ratio  (<  1) 

3.  Effective  reflectance  of  target  (<  1) 

4.  Effective  emittance  of  target  (<  1) 

5.  Total  radiant  intensity  (using  Fresnel  equations) 

6.  Total  radiant  intensity  (assuming  uniformity) 

7.  Aspect  angle  (degrees) 

8.  Error  in  using  normal  reflectance 

9.  Error  in  using  normal  emittance 

10.  Nadir  angle  (degrees) 
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TABLE  4 

INPUT  DATA  FOR  DIREC 


TJ 

QJ  U 

03  Q) 

ZD  CO 

M-l 

<U  -H  4J 
P  O 

1  z 

> 


180 


TABLE  4  (Cont.) 
INPUT  DATA  FOR  DIREC 
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TABLE  4  (Cont.) 
INPUT  DATA  FOR  DIREC 
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VI1*  CQDE  FQR  POLARIZATION  CONTENT  OF  TARGET  EMISSION  SIGNATURE  (EP1) 

A.  INTRODUCTION 

Computer  code  EP1  was  designed  1  calculate  the  polarization  con¬ 
tent  of  the  emission  component  of  a  target's  optical  signature.  Input 
tables  from  the  computer  code  TARG  (Appendix  I)  are  required,  so  assump¬ 
tions  implicit  in  it  are  pertinent  here.  The  general  procedure  is  to 
(1)  subdivide  the  target  into  elemental  areas,  (2)  determine  the  two 
components  of  polarized  emissivity  for  a  general  element,  (3)  reference 
these  components  to  a  plane  independent  of  the  particular  element,  and 
(4)  integrate  over  the  portion  of  the  target  exposed  to  the  detector. 

It  is  assumed  that  these  elements  may  be  considered  flat  planes  and  that 

the  target  is  small  relative  to  its  range  from  the  detector. 

B.  DEVELOPMENT  OF  EQUATIONS 

If  the  latitude  \p,  longitude  A,  and  altitude  3  (normalized  to  the 
radius  of  the  earth)  of  the  target  (denoted  by  subscript  t)  and  of  the 
detector  (subscript  d)  are  given,  their  position  vectors  R  may  be  written 
with  respect  to  the  earth-fixed  coordinate  system  shown  in  Fig.  13,  as: 

a  A 

Rt  =  ^t^ex  C°S  ^t  COS  Xt  +*6y  COS  ^t  Sin  *t  +  ez  Sin  ^  (VII-l) 

~y  .a  A 

Rd  =  3d^ex  COS  ^d  C0S  Xd  +  ey  C0S  ^d  Sin  Xd  +  &z  Sin  ^  (VII_2) 

and  the  vector  from  the  target  to  the  detector  is  simply: 

r  =  Rt  -  Rd  (VI 1-3) 

and 

r  "  r/|r|  (VII-4) 
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A 


Figure  13.  Earth-Fixed  Coordinate  System 
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Alternatively,  if  the  nadir  angle  6  and  the  azimuth  angle  x  to  the 
viewer  relative  to  the  target  are  specified  (see  Fig.  14)  we  can  write 
simply: 


r  -  e  sin  6  cos  x  +  e  sin  6  sin  x  “  e  cos  6 
*  y  z 


(VII-5) 


where  the  target  is  located  on  the  z  axis. 

Consider  an  arbitrary  element  of  the  target  whose  size  and  orien¬ 
tation  is  specified  by  a  normal  vector  A.  The  polarized  components  of 
radiant  intensity  from  this  element  are: 


dIp  =  ^nbb  dA  cos  9 


dIL  =  lNBB  dA  CoS  9 


(VII-6) 


where  ep  and  eL  are  the  parallel  and  perpendicular  polarization  sensitive 
emissivities,  respectively,  N^B  is  the  blackbody  radiance  of  the  element, 
dA  is  the  area  of  the  element  and  9  is  the  angle  between  the  element's 
normal  vector  r.  Using  vector  algebra,  the  angle  0  is  given  by: 

cos  0  =  r  •  A  (VII-7) 

Using  this  angle,  theoretical  polarized  emissivities  e  and  e  are  calcu- 

P  L 

lated  by  Fresnel  equations.  This  polarization  is  defined  with  respect  to 
the  plane  of  incidence  of  the  particular  target  element  being  considered 
and  so  must  be  referenced  to  a  common  plane.  If  tp  is  the  angle  between 
the  plane  of  incidence  (whose  unit  normal  is  M)  and  an  arbitrary  reference 
plane  containing  r  (whose  unit  normal  is  P) ,  the  polarization  components 
of  radiant  intensity  with  respect  to  the  equatorial  plane  are  simply: 
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dip  =  dip  cos^  4>  +  dl^  sin^  <p 
dl^  =  dip  sin^  <J>  +  dl^  cos^  <p 


(VII-8) 


Note  that  since  the  sine  and  cosine  functions  are  squared,  it  does  not 
matter  whether  the  angle  <p  or  its  supplement  (it  -  cp )  is  used.  From  vec¬ 
tor  algebra 


cos  <f>  =  M  •  P  (VII-9) 

*f 

where  M  may  be  written  as: 

M  =  A  *  r  and  M  -  M/|m|  (VII-10) 

As  an  additional  check  on  the  calculation,  an  option  is  included  to  allow 
selection  of  any  reference  plane  (which  includes  r) .  The  control  param¬ 
eter  is  NREF  where,  if  NREF  «  1,  the  plane  is  defined  by  vectors  r  and 
R. ,  or,  if  NREF  >  1,  the  vectors  r  and  Q  are  used  (Q  is  an  arbitrary 
vector  to  be  set).  Thus  the  reference  plane's  normal  vector  is: 

P  =  Rt  x  r  if  NREF  £  1  (VII-11) 


or 


P  =  Q  X  r  if  NREF  >  1  (VII-12) 

For  computer  evaluation,  Eqs.  VII-6  with  Eq.  VII-8  are  converted  to 
finite  difference  form,  numerically  integrated,  and  normalized  with  re¬ 
spect  to  N„t  and  A  : 

BB  p 
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4  cos  ^ 

i 

,  2 
(  £p  cos 

*  + 

.  2 
i  sin 

M/Ap 

(VII-13) 

4  =£aa  c°s  0, 

i  1 

(  £p  sin 

<j>  + 

2 

i  COS 

«/Ap 

(VII-14) 

Finally,  these  results  are  expressed  in  terms  of  a  parameter  commonly 
used  as  a  measure  of  the  polarization: 


e 


P 


i  -  i 
i+  i 


(VII-15) 


It  is  seen  that  £p  is  independent  of  the  target  temperature  if  the  target 
is  assumed  isothermal,  (i.e.,  N  uniform). 


As  a  simple  check  on  the  code,  the  projected  area  A  of  the  target 

P 

is  also  calculated.  This  may  be  written  as: 


A 

P 


cos  0  dA 


and 


A  AA.  cos  0. 

P  ^  i  i 


(VII-16) 


(VII-17) 


where  the  only  elements  considered  are  those  for  which  the  cosine  terms 
are  positive. 

C.  INPUT/OUTPUT  SPECIFICATIONS 

Before  EP1  can  be  called,  TARG  must  be  requested  by  the  user  (see 
Appendix  I) .  As  many  calls  as  required  may  then  be  made  to  EP1  without 
recalling  the  TARG  code.  Table  5  defines  the  input  necessary  for  using 
this  code.  The  first  data  card  provided  for  input  to  EP1  must  contain 
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the  letters  EP1  in  columns  one  through  three  to  indicate  the  operator 
desires  to  use  this  code.  The  output  generated  by  EP1  follows: 

1.  Aspect  angle  (degrees) 

2.  Nadir  angle  (degrees) 

3.  Azimuth  angle  (degrees) 

4.  Projected  area  (sq  mi) 

5.  Degree  polarized 

6.  Parallel  polarization 

7.  Perpendicular  polarization 
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TABLE  5 

INPUT  DATA  FOR  EP1 
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TABLE  5  (Cont.) 
INPUT  DATA  FOR  EP1 


TABLE  5  (Cont.) 
INPUT  DATA  FOR  EP1 
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VIII.  CODE  FOR  POLARIZATION  CONTENT  OF  SPECULARLY  REFLECTED  EARTHSHINE  iSPOL) 

A.  INTRODUCTION 

This  computer  code  is  designed  to  calculate  the  polarization  content 
of  the  specularly  reflected  component  of  a  target's  optical  signature. 
Theoretical  Fresnel  expressions  are  used  to  provide  the  polarized  reflec¬ 
tivities.  Input  tables  calculated  by  TARG  are  required,  so  the  assumptions 
implicit  in  this  code  are  also  pertinent  here.  The  general  procedure  is 
to  (1)  subdivide  the  target  into  elemental  areas,  (2)  determine  the  two 
components  of  polarized  reflectivity  for  each  element,  (3)  reference 
these  components  to  a  plane  independent  of  the  particular  element,  and 
(4)  integrate  over  the  portion  of  the  target  exposed  to  the  detector  that 
will  specularly  reflect  earthshine.  It  is  assumed  that  these  elements 
may  be  considered  flat  plates  and  that  the  target  is  small  relative  to 
its  range  from  the  detector. 

B.  DEVELOPMENT  OF  EQUATIONS 

The  polarized  components  of  radiant  intensity  of  an  incremental 
area  dA  of  the  target  specularly  reflecting  radiance  L.^ C 4>)  from  an  earth 
element,  toward  a  detector  (Fig.  3)  can  be  written: 

dl  =  P?L,($)  dA  cos  0  (VIII-1) 

L  LA 

dl  =  p^L.  (<f>)  dA  cos  0  (VIII-2) 

r  FA 

s  s 

where  pp  and  p^  are  the  parallel  and  perpendicular  polarization  sensitive 
reflectivities,  respectively. 

To  obtain  the  total  polarized  radiant  intensity  due  to  reflected 
earthshine,  this  expression  is  integrated  over  all  the  target  that  is 
visible  to  the  detector  and  from  which  specular  reflection  from  the  earth 
is  possible. 
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The  angle  0  Is  defined  in  Fig.  3  and  is  determined  as  follows: 

R  =  V  -  T  (VIII-3) 

cos  0  ■  R  •  n/Rl  (VIII-4) 

where  n  is  the  normal  to  the  target  element  and  R1  is  the  magnitude  of 

->• 

R  and  is  obtained  from 

R!  -  (Rx  +  +  Rz)1/2  (VIII-5) 

where  R  ,  R  ,  and  R  are  the  components  of  R  in  an  earth  fixed  coordinate 
x  y  z 

system. 

If  0  is  greater  than  90°,  there  is  no  radiation  to  the  viewer. 

When  0  is  less  than  90°,  one  must  for  each  element  determine  the  vector 
i  and  ensure  that  it  intersects  the  earth.  Three  independent  equations 
are  required  in  order  to  determine  the  three  components  of  i.  Two  inde¬ 
pendent  equations  are  determined  from  the  relationship  (Fig.  3) 

h  x  t  =  £  x  ft  (VIII-6) 

and  the  third  from  the  relationship 

f  •  n  =  n  •  i  (VIII-7) 

To  ensure  that  the  vector  i  intersects  the  earth,  the  angle  6 
(Fig.  3)  is  calculated  from  the  relationship 
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T  •  i  =  cos  0 


(VIII-8) 
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If  3  <  4>,  where  <j)  =  arc  sin  (earth  radius/earth  radius  and  target 
altitude),  then  the  vector  i  intersects  the  earth  and  the  element  contri¬ 
butes  to  the  specularly  reflected  radiant  intensity  of  the  target. 

The  polarization-sensitive  reflectivities  are  determined  from 
Fresnel's  equations  with  respect  to  the  plane  of  incidence  of  the  par¬ 
ticular  target  element  being  considered  and  must  be  referenced  to  a 
common  plane. 

The  common  plane  is  the  plane  determined  by  the  vector  from  the 
center  of  the  earth  to  the  target  and  the  vector  from  the  center  of  the 
earth  to  the  detector.  If  these  vectors  are  collinear,  the  common  plane 
is  determined  by  these  vectors  and  the  vector  from  the  center  of  the 
earth  through  the  North  Pole.  If  all  three  of  these  vectors  are  collinear, 
the  common  plane  is  the  plane  determined  by  these  vectors  and  a  vector 
from  the  center  of  the  earth  through  the  equator  at  the  Greenwich 
Meridian.  The  polarization  components  of  radiant  intensity  with  respect 
to  the  common  plane  is  simply 


dIL 

K) 

2 

|  COS  <j)  + 

(dIp) 

•  2  A 

sin  (p 

(VIII-9) 

dIp  = 

K) 

2 

sin  0  + 

K) 

2  A 

cos  cp 

(VIII-10) 

where  <j>  is  the  angle  between  the  plane  of  incident  and  reflected  radiation 
and  the  common  plane.  Note  that  since  the  sine  and  cosine  functions  are 
squared,  it  does  not  matter  whether  the  angle  cp  or  its  supplement 
(n  -  <f>)  is  used.  From  vector  algebra 
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cos  <J>  =  m  •  p 


(VIII-11) 
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where 


i  s  f  x  n 


(VIII-12) 


and  p  is  the  unit  normal  vector  to  the  reference  plane. 

P  =  i  x  Q  (VIII-13) 

where  Q  is  one  of  the  three  vectors  indicated  earlier. 


For  computer  evaluation,  Eqs.  VIII-9  and  VIII-10  with  Eqs.  VIII-1 

and  VIII-2  are  converted  to  finite  difference  form,  numerically  integrated, 

and  normalized  with  respect  to  L. ($)  and  A  : 

A  p 


V 


COS 


2 

cos 


<fi  +  Pp  sin2 


(VIII-14) 


:p  =X>i  cos 


S  .  2 
PL  sin 


,  s  2 
+  pp  cos 


(VIII-15) 


Finally  these  results  are  expressed  in  terms  of  a  parameter  commonly 
used  as  a  measure  of  the  polarization: 


I 


f 

I 

L 


» 


t 


(VIII-16) 


As  a  check,  the  projected  area  (A  )  of  the  target  is  calculated. 

It  is  simply 

Ap  =EAAi  Cos  e±  (VIII-17) 

where  the  only  elements  considered  are  those  for  which  the  cosine  terms 
are  positive. 
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C.  INPUT/OUTPUT  SPECIFICATIONS 

Before  SPOL  can  be  called,  TARG  must  be  requested  by  the  user  (see 
Appendix  I).  As  many  calls  as  required  may  then  be  made  to  SPOL  without 
recalling  the  TARG  code.  Table  6  defines  the  input  necessary  for  using 
this  code.  The  first  data  card  provided  for  input  to  SPOL  must  contain 
the  letters  SPO  in  the  first  three  columns  to  indicate  the  operator 
desires  to  use  this  code.  The  output  generated  by  SPOL  follows: 

1.  Aspect  angle  (degrees) 

2.  Nadir  angle  (degrees) 

3.  Azimuth  angle  (degrees) 

4.  Projected  area  (sq  mi) 

5.  Degree  polarized 

6.  Parallel  polarization 

7.  Perpendicular  polarization 
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INPUT  DATA  FOR  SPOL 
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TABLE  6  (Cont.) 
INPUT  DATA  FOR  SPOL 
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TABLE  6  (Cont.) 
INPUT  DATA  FOR  SPOL 
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o 

o 

00 

rH 

O 

O 

00 

O 

o 

00 

rH 

• 

0 

o 

o 

00 

00 

rH 

rH 

O 

O 

4J 

4J 

• 

. 

O 

O 

CO 

CO 

44 

0) 

44 

0) 

u 

S  X 

00 

4  00 

QJ 

41 

XJ 

T3 

eg  4j 

oj  u  tn 

C  X 
O  44  -H 

X  4-1 

c 

eg  O  44 

X  x  x 


eg  4-i 
44  o  tn 
x  eg 

Q  HJ  H 

8  xi 

4J 

4-1  o  X 


^  iH 
«  0 

H  X 

M 

&  B 

Q 

Q 

P  ^ 

fcw  m 

C  cn 

O  44  tg 

X  H 

e 

eg  o  tu 
XXX 
X  X 

tn 

44  -H  tn 
X  *H 

g  4i 
6  H  tn 
60  -H 

4-t  C  X 

H  eg  x 
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APPENDIX  I 

TARGET -ELEMENT  CODE  (TARG) 


A.  INTRODUCTION 

This  program  is  one  of  a  series  to  be  employed  for  calculating  the 
radiant  intensity  of  exoatmospheric  targets.  Its  specific  purpose  is  to 
determine  the  area  normal  vectors  of  elemental  areas  on  the  target.  The 
components  of  these  vectors  are  stored  for  use  in  subsequent  codes.  Only 
sphere-cone-shaped  targets  may  be  handled  by  this  program.  However,  the 
dimensions  and  cone  angles  are  adjustable  so  the  shapes  can  include  pure 
cones,  hemispheres,  and  hemisphere-cylinders.  The  base  of  these  targets 
is  assumed  to  be  a  simple  circular  plate. 


Position  of  the  target  must  be  supplied  in  terms  of  its  latitude 
and  longitude.  The  orientation  of  the  target  is  specified  by  two  other 
angles.  Components  of  the  required  vectors  are  given  in  the  earth-fixed 
coordinate  system  defined  by  the  North  Pole  and  Greenwich  meridian. 

B.  DEVELOPMENT  OF  EQUATIONS 

An  earth-fixed  Cartesian  coordinate  system  (x,  y,  z)  is  defined 
by  the  origin  at  the  center  of  the  earth,  the  z-axis  pointing  north  and 
the  x-axis  toward  the  Greenwich  Meridian  (all  coordinate  systems  are 
right-handed) .  Components  of  the  target  position  unit  vector  T  T  T 

X  *  y  *  2 

(see  Fig.  15)  are  given  in  terms  of  its  longitude  XT  and  its  latitude 
(positive  in  Northern  Hemisphere,  negative  in  Southern),  by: 

Tx  =  cos  cos  XT 

Ty  =  -cos  sin  XT  (AI-1) 

T  =  sin  \p 
z  I 
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A  local  coordinate  system  (x^  z J  ±s  defined  such  that  the 

Xi  3X18  18  colllnear  with  Tx  and  zQ  axis  in  the  x^  plane  (see  Fig.  16a). 
This  transformation  may  be  expressed  by  the  following  equations  where 

Xl’  yl’  Z1  are  components  of  any  vector  in  the  local  coordinate  system 

and  xq»  yQ»  zq  are  components  of  the  same  vector  in  the  earth-fixed 
system: 

*  =  xl  cos  COS  XT  ~  yi  sin  ~  zi  siu  ^x  C0S  XT 

y  =  xi  cos  sin  XT  +  yi  COS  Xx  ~  Z1  sin  ^x  Sin  XT  (AI-2) 

z  =  x^  sin  ii»T  +  z^  cos  ip 

A  body-fixed  coordinate  system  y^  z2  is  defined  auch  that  the 
X2  axis  lies  along  the  vehicle  axis  from  its  base  toward  its  nose.  The 

z2  axis  is  set  arbitrarily  to  some  direction  consistent  with  the  value 

of  the  material  orientation  angle  used.  Vehicle  orientation  with  respect 
to  the  local  coordinate  system  is  given  by  two  angles  0T  and 
Transformation  is  accomplished  (see  Fig.  16b)  by  first  rotating  the  x 
yl»  zx  system  about  the  y^  axis  (counterclockwise  when  looking  along  y  ) 
an  amount  ^  so  that  moves  to  z^  This  system  is  then  rotated  about 
the  z2  axis  (clockwise  when  looking  along  z^)  an  amount  0T  so  that  x 
and  yx  move  to  x2  and  y2  respectively.  This  transformation  may  be 
expressed  by  the  following  equations,  where  x2>  y£ ,  z£  are  components  of 
any  vector  in  the  body-fixed  system  and  x^  y±,  Z]^  are  components  of  the 
same  vector  in  the  local  system: 

X1  *  X2  COS  0T  C0S  ^T  ~  y2  sin  9t  C0S  ^T  +  z2  Sin  ^T 

yX  =  x2  Sln  9t  +  y2  COS  6T 

Z1  =  ~X2  COS  0T  Sln  ^T  +  y2  Sin  0T  Sin  ^T  +  z2  COS  *T 
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(AI-3) 
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Combining  Eqs.  AI-2  and  AI-3,  it  is  possible  to  express  the  components 

x,  y,  z  in  terms  of  x  .  v  .  z  : 

2  3  2  2 

x  =  x2(cos  0T  cos  cos  cos  XT  -  sin  eT  sin  AT 

+  cos  0  sin  <|>  sin  ip  cos  A)  -  y  (sin  0^  cos  d>  cos  iL 

z  T  T  T 

cos  Ax  +  cos  0T  Sin  AT  +  sin  0T  sin  ^  sin  ip^,  cos  A  ) 

+  z^  (sin  <j>r£  cos  ipT  cos  AT  -  cos  4>T  sin  ipT  cos  A^)  (AI-4a) 

y  =  x2(cos  0T  cos  4>t  cos  PT  sin  AT  +  sin  0T  cos  AT 
+  cos  0T  sin  4>t  sin  ^  sin  XT)  +  y2(cos  0T  cos  A 

-  sin  0T  cos  <^T  cos  iPT  sin  A?  -  sin  @T  sin  <J>T  sin  ^ 

sin  At)  +  z2(sin  4>T  cos  ^  sin  AT  -  cos  d>T  sin  ^  sin  AT 

,  (AI-4b) 

z  x2(cos  @x  cos  <j>x  sin  ^  -  cos  0^  sin  <p^  cos  ip  ) 

+  y2(sin  0T  sin  <f>T  cos  ^  -  sin  0T  cos  <J>T  sin  ip  ) 

+  z2(sin  <pT  sin  ^  +  cos  ^  cos  tf^)  AI-4c) 

Elemental  areas  on  the  targets  are  defined  (see  Fig.  17)  by  sub¬ 
dividing  the  spherical  portion  angularly  in  two  directions,  <f,  and 

9S*  the  cone  an8ularly  in  one  direction,  6^  and  the  base  as  one  ele¬ 
ment.  Components  of  these  vectors  in  the  body-fixed  coordinate  system 
are  given  by: 

89 

UNCLASSIFIED 


asphere 


f 
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t 


\  =  \  =  0 
X  y 


\  =0 


(AI-5) 


Ar  “  Ar-  sin  a  K  =  K  cos  a  sin  6  A  =  A  cos  a  cos  6  (AI-6) 

^  L  c  c  CCC  c  '  / 

y  z 


X 


AS  “  AS  C0S  ^S  AS  =  AS  sin  *s  sin  0-  A-  =  A-  sin  cos  e-  (AI“7) 


s  s  s 
z 


where  A^  - 


A  =  A0  A<f)  R  sin  <p 
s  s  s  n  s 


(AI-8) 


A  = 
c 


A0 

c  ,..2  .2  2  , 

~n  a  (Rb  "  Kn  C°S  a) 


or  if  a  -  0,  A  -  A0  R  L,  where  a  is  the  cone  half  angle,  R  is  the  nose 

C  C  II  n 

radius,  R^  is  the  base  radius,  and  L  is  the  cylinder  length.  An  integer 

number  of  angular  elements  are  specified  (N  ,  N  ,  N  )  so  that: 

(pud) 
s  s  s 


46 

c  Ne 


46 

S  N8 


A<PC  =  (j  ~  a)/N 

S  £.  q 


(AI-9) 


As  a  check  on  the  program,  the  total  projected  area  of  the  target 
may  be  calculated  given  the  observer  unit  position  vector  components 
^x»  ^y»  the  earth-fixed  system.  This  vector  is  dotted  with  each 

area  normal  vector  to  give  the  projected  area  that  each  element  presents 
to  the  viewer: 


A  =  A  cos  v  =  A  •  V 
P 


=  A  V  +  AV  +  AV 
xx  y  y  ?  z 
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If  the  dot  product  is  negative,  indicating  the  area  if  facing  away  from 
the  viewer,  the  element  is  not  considered.  All  the  other  contributions 
are  summed  to  give  the  total  projected  area  of  the  target. 


C.  INPUT/OUTPUT  SPECIFICATIONS 

Table  7  defines  the  input  necessary  to  use  this  code.  The  first 
data  card  must  contain  the  letters  TAR  to  indicate  the  operator  desires 
to  use  this  code.  All  codes  except  SURAD  require  the  output  generated 
by  TARG. 


Special  target  shapes  may  be  handled  by  proper  selection  of  R^ ,  R^, 

L  and  a.  A  pure  cone  is  specified  by  setting  R^  to  zero,  a  hemisphere 

by  setting  a  and  L  to  zero,  a  hemisphere-cylinder  by  setting  a  to  zero 

and  simple  circular  flat  plate  (one  sided)  by  setting  R  ,  a,  and  L  to 

n 

and  Nrt  should  also  be  set  to  zero,  and 


zero. 


If  R  is  set  to  zero,  N, 

n  <fi 


0 


similarly  for  NQ  if  XL  is  zero. 

0 


Results  (the  X-array)  are  written  on  a  disc  file  labeled  TAPE4  for 
use  in  latter  programs.  Aside  from  intermediate  output  generated  by 
TARG  for  debug  purposes  when  W  equals  "1."  or  "2.",  no  output  is  generated 
for  the  printer. 


*X(I) 


X(j  +  1),  X(j  +  2),  X(j  +  3)  = 

A  ,  A  ,  A 
x  y  z 

X(k  +  1) ,  X(k  +  2) ,  X(k  +  3)  = 

M  ,  M  ,  M 
x  y  z 

where 


j  =  1,  2,  ...,  N 
k=N+l,  N  +  2,  N  +  3 
. . . ,  2N 

N  =  NPS  •  NTS  +  NTC  +  1 


*OUTPUT 
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TABLE  7 

INPUT  DATA  FOR  TAR6 


TABLE  7  (Cont.) 
INPUT  DATA  FOR  TARG 
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o 

II 


« 
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Both  NPS  and  NTS  equal  zero  if  RN 


nge  of  I Value  Used 


UNCLASSIFIED 


APPENDIX  II 

EARTH-ELEMENT  CODE  (EARTH) 


A.  INTRODUCTION 

This  digital  computer  code  is  one  of  a  series  designed  to  calculate 
the  radiant  intensity  of  exoatmospheric  objects.  The  specific  purpose 
of  this  code  is  to:  (1)  subdivide  the  portion  of  the  earth  which  irradi¬ 
ates  a  given  target  location  into  a  specified  number  of  elemental  areas, 
(2)  calculate  the  components  of  the  unit  vectors  pointing  toward  this 
target  position  from  each  elemental  area  on  the  earth,  and  (3)  calculate 
the  irradiance  at  the  target  location  due  to  b)th  earth-emitted  radiance 
and  earth-reflected  sunlight  given  the  posit  ion  of  the  sun.  These  results 
are  used  in  subsequent  computer  codes,  DIFUS1  and  TINTF. 

B.  DEVELOPMENT  OF  EQUATIONS 

Components  of  the  vectors  are  found  in  an  earth-fixed  coordinate 
system  (x,  y,  z)  defined  (see  Fig.  15)  by  the  origin  at  the  center  of 
the  earth,  the  z-axis  pointing  north  and  the  x-axis  toward  the  Greenwich 
Meridian  (all  coordinate  systems  herein  are  right-handed  Cartesian).  The 
target  location  vector  T  in  this  earth-fixed  coordinate  system  can  be 
given  in  terms  of  the  target's  longitude,  AT,  its  latitude,  ifj  (positive 
in  the  Northern  Hemisphere,  negative  in  the  Southern  Hemisphere),  and 
its  altitude  0  (normalized  to  the  earth's  radius),  by: 


T  =  iS  cos  ipT  cos  At  -  j0  cos  i|i  sin 


XT  +  k0  sin 


(AII-1) 


where  i,  j,  and  k  are  unit  vectors  along  the  x,  y,  and  z  axes  respec¬ 
tively.  Using  the  sun's  longitude  A  and  latitude  \ p  ,  the  unit  vector 

s  s 

toward  the  sun  is  similarly: 
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S  “  i  cos  Ip  cos  X  -  j  cos  if  sin  X  +  k  sin  \p  (AII-2) 

s  s  s  s  s 

For  convenience  a  local  coordinate  system  (x.^ ,  y^ ,  z )  is  set  up 
with  the  origin  again  at  the  center  of  the  earth,  the  z^  axis  pointing 
toward  the  target  and  the  z  axis  in  the  ^1~z1  plane.  This  transforma¬ 
tion  may  be  expressed  by  the  following  equations,  where  x^ ,  y  ,  z^  are 
components  of  any  vector  in  the  local  coordinate  system  and  x,  y,  z 
are  components  of  the  same  vector  in  the  earth-fixed  system: 

x  =  x^  sin  cos  X^,  +  y^  sin  X^  +  z^  cos  cos 
y  *  -x1  sin  ij»T  sin  XT  +  y1  cos  XT  -  z±  cos  \p  sin  XT  (AII-3) 
z  =  — X1  cos  ipT  +  z 1  sin  1/1 

Assuming  a  spherical  earth,  the  portion  which  irradiates  the  target 
(see  Fig.  18)  is: 

cos  =  1/8  (AII-4) 

This  earth  disk  is  subdivided  into  n^  elemental  angles  in  the  <j>  direction 
and  n^  in  the  0  direction,  so  that  the  sizes  of  the  elements  are 

A<p  =  <p  / n 
m  ® 

A0  =  2Tr/nQ  (AII-5) 

AA  =  A<J>A0  sin  <p  , 
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The  area  normal  vector  A,  having  a  magnitude  equal  to  the  size  of  the 
element  and  a  direction  normal  to  it,  may  be  expressed  in  the  local 
coordinate  system  by: 

A  =  (i  sin  <p  cos  0  +  j  sin  $  sin  6  +  k  cos  $)  A*A0  sin  <p  (AII-6) 

The  position  of  each  of  the  elemental  areas  is  given  by  a  unit  vector 
E,  expressed  in  the  local  system  by: 


E  -  i  sin  <f>  cos  0  +  j  sin  <f>  sin  0  +  k  cos  <p  (A 

Both  these  vectors  are  transformed  to  the  earth-fixed  coordinate  system 
using  Eq.  AII-3.  By  vectorially  subtracting  E  from  T  we  get  the  range 
vector  R  from  the  area  element  to  the  target: 


(AII-8) 


Assuming  the  elemental  areas  radiate  diffusely,  the  irradiance  at 
the  target  due  to  any  one  of  them  is: 


EE  *  (NA/R2)  cos  yet 


(AII-9) 


where  N  is  the  radiance  of  the  earth  element,  yet  is  the  angle  between 
A  and  R.  Spectral  earth  radiance  at  any  wavelength  is  assumed  dependent 

only  on  the  local  nadir  angle  "n"  (see  Fig.  19),  which  can  be  expressed 
as : 


n 


tan 


-1 


sin  j> 

3  -  cos  <p 


(AII-10) 
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I 


CENTER  OF 
EARTH 


Figure  19.  Sun-Earth-Target  Geometry 
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M  =  s  x  A  and  N  =  A  x  r  (AII-15) 

An  option  is  included  as  a  check  nn  the  code  to  calculate  the  total 
view  factor  between  the  target  and  the  portion  of  the  earth  irradiating 
it,  for  emitted  earthshine  Fp  and  for  reflected  sunlight  F  . 

C.  INPUT/OUTPUT  SPECIFICATIONS 

Tables  8  and  9  define  the  input  necessary  to  use  this  code. 

Table  9  is  needed  only  if  earth-reflected  sunlight  is  to  be  considered. 
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The  first  data  card  provided  after  the  TARG  input  data  must  contain  the 
letters  EAR  to  indicate  the  operator  desires  to  use  this  code.  Only  codes 
TINTF  and  DIFUS1  require  the  output  from  EARTH. 

Output  is  written  in  three  sets.  The  first  and  second  write  the 
x,  y,  and  z  components  of  the  R  vector,  E_  and  E  for  each  earth  element. 

L  S 

However,  to  save  storage  and  writing  time,  if  the  sun  does  not  irradiate 

an  element,  the  sign  of  E_  is  reversed  and  no  value  of  E  is  written.  The 

fc.  s 

first  set  is  written  in  binary  on  a  disc  file  labeled  TAPE10  for  use  in 
subsequent  computer  codes  TINTF  and  DIFUS1,  and  the  second  (if  WRITE  is 

set  equal  to  one)  is  printed  out  for  debug  purposes.  In  these  outputs, 

the  x,  y,  and  z  components  (in  that  order)  of  all  the  earth  elements  are 

written  first  in  the  same  order  in  which  they  were  calculated;  that  is, 

all  the  0's  in  the  first  <p  ring,  then  all  the  0's  in  the  second  <f>  ring, 

and  so  on,  until  all  the  4>  rings  are  written  After  the  components  of  the 

R's,  the  values  of  E  and  E  are  written  in  the  same  order  as  above. 

ii  s 

First  the  E  then  the  E  for  each  earth  element  (unless  E„  is  negative, 
iii  s  hi 

in  which  case  Eg  is  not  given).  The  third  set,  the  "point  view  factors" 
which  the  total  and  the  sunlit  portions  of  the  earth  present  to  the  target, 
are  printed  out  if  CK  is  set  equal  to  one. 
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APPENDIX  III 

EMITTED  RADIATION  CODE  (HRAD) 


A.  INTRODUCTION 

This  code  is  designed  to  calculate  the  radiant  emission  of  a  target 
given  to  temperature  and  emissivity.  The  radiation  emitted  by  a  reentry 
vehicle  is  a  function  of  its  surface  temperature  and  spectral  emissivity. 
For  the  purpose  of  numerical  calculation,  the  surface  is  divided  into  "n" 
sections  which  are  sufficiently  small  to  be  assumed  isothermal.  For  a 
Lambertian  radiator,  the  radiant  intensity  is 


C  n=N 

■tE 

0 


A  e(A,  T  )  A 
n  n 


-5[» 


(C./AT  ) 


where 


2 

A  is  the  projected  area  (cm  ) 

1(\\  is  the  radiant  intensity  (W-sr_1) 


e(A,  T)  is  the  spectral  emissivity 
A  is  the  wavelength  (pm) 

T  is  the  surface  temperature  (°R) 
C1  =  3.7413  x  104  W/cm2-um4 

C2  =  2.59  x  io4  pm-0!* 


The  radiant  intensity  is  obtained  by  integrating  over  wavelength: 


C.  n=N  r 

■tZ >J 

r\  •'  i 


AA 


c(A, 


T  )  A 
n 


-5 


(C,/AT  ) 
l  n 


dA 
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A  series  approximation  to  this  integral  of  Planck's  equation  that  is  used 
in  the  program  follows: 


where  X.  =  C./A.T 
1  2  1  n 


C  /A  T 
2  2  n 


The  HRAD  program  calculates  the  radiance  and  radiant  intensity, 
spectral  and/or  integrated,  for  the  surface  in  question,  given:  (1)  tem¬ 
perature,  emissivity,  and  size  of  isothermal  areas;  (2)  spectral  wave¬ 
lengths;  and  (3)  limits  of  integration. 

The  assumptions  incorporated  in  the  program  are: 

1.  Multiple-isothermal-area  concept 

2.  Lambertian  radiator 

3.  Emissivity  a  function  of  area  and  wavelength 

The  latter  assumption  does  not  reflect  the  difficulty  inherent  in 
the  inclusion  of  variable  emissivity  but  rather  a  lack  of  information  on 
the  spectral  and  thermal  behavior  of  optical  properties  of  current  reentry 
vehicle  heat  shield  materials. 
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B.  COMPUTER  CODE 


T  HF  HR  A  n  PROGRAM  CALCULATES  THE' 
SRFCTRAL  AMD/QR  INTEGRATED.  FOR 
1.  TFYPFPA  tijrf  ,  fmisstvity. 

SPFPTRAL  WAVELENGTHS 
A.  LIMTTS  Of  INTEGRATION 
TRF  ASSUMPTIONS  INCORPORATED  IN 
1.  NIJLT1PLE-ISOTMFRMAL-AREA 
P.  LAMDrRT  IAN  RADIATOR 

P.  RMISSIVI^Y  a  function  OF 


RADIANCE  AND  RADIANT  INTFNSITY. 
THF  SURFACF  IN  QUESTION.  GIVEN 
AND  SIZE  OF  ISOTHERMAL  AREAS 


THE  PROGRAM  ARF 
CONCEPT 

ARFA  AND  INTEGRATION  WAVELENGTH 


TNfut  SETUd 

CARO  1  -  IDENTIFICATION  TO  HF  PRINTED  ON  THF  RIJN  (FORMAT  1 6 AS > 
CAPO  ?  -  N  A  SPEC .  NARFAS.  NT  IMPS.  NLAMRS.  NOLAMS  (FORMAT  5110) 

NAGDfC  -  numhfr  op  ASPECTS 
NARFAS  -  NIIMRFP  OF  APF  A  S 
NTTMFS  -  NIIMRFR  OF  TIMES 

M  L  A  m  q  c;  -  NIIMHFR  QF  SPfcTPAL.  WAVFLENGTHS 
NOLAMS  -  NIIMHFR  OR  INTEGRATED  WAVELENGTH  RANDS 


CARO  SF  T  0 


(  A G R E o ( T ) , 1  =  1 .NASPFC) 


ardfc 


THE  valiifs  OF  THE  ASPECT  ANGLES 


FORMAT  (HF10.6) 


CA&n  SET  4 
NT  YDE  THF 


DESCRIPTION  OF  ARFAS.  THE  INFORMATION  WILL  VARY 

rrpr  rus?tJrFJHF'  TVPE  °F  SURFacE  “ITH™  section 


N  T  YDE  = 1 
I N  a  U  T  IS 
RAD 
ANSI  F 


SPHERICAL  CAR 

UTYPF,  pan,  ANGLF.  FORMAT  (I5.PE10.3) 

GAP  RADIUS 

THF  ANGLE  FROM  the:  CENTFR  OF  THF  HEMISPHERE  TO  THE  rnrr 
OS  THE  C»=.  IF  THF  CP  IS  A  HEMISPHERE.  ANGLcJ 


N  T  Y  D  F  -?  CYLINDER 

I\nuT  I S  NTYPF,  RAD.  LENGTH  FORMAT 
RAn  -  THF  CYLINDER  RADIUS 

Length  -  the  length  of  thf  cylinder 


(  IS.PE10. 3) 
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NTYPF  =3  0®  4  CONF  OP  CONIC  FUNCTION 
INPUT  IS 

1.  NTYPF. *  PE  AR  RADIUS  OF  VFH  I CLF  ♦  NO.  OF  CONFS*  CONE  ANGLE 
FORMAT  (TS.3F10.3) 

?.  XLFNU-NO.  OF  CONFS)  FORMAT  ( BF  10.6) 

Xlfm  -  JHF  LENGTHS  FOd  each  CONE  FRUSTRUM 

NTYPF  =5  FLAT  PLATE 

Input  is  NTYPE*  RAO,  ANGLF  FORMAT  (15* ?F 10.3) 
pftn  -  THF  RADTIJS  OF  THE  PLATE 

ANGLE  -  THF  ANGLF  RETWFEN  THE  PLATE  NORMAL  AMD  THE  AXIS  OF 
SYMMETRY.  USUALLY  WILL  RE  FITHF»  0  OR  1 A 0  DEGREES 

NTYDF  =5  REAR  SDHFRICAL  PLATF  FOR  INPUT  SFF  NTYPE=1 


CAPO  SFT  5  -  FMISS('IDI  AMS*NARFAS)  (FORMAT  BF10.6) 

FMTTTANCS  FOR  EACH  area  AND  FOP  EACH  INTFGRATFD  WAVELENGTH. 
EACH  OARD  RFPRFSFNTF  PMTTTANCE  VALUFS  FOR  A  WAVELFNGTH  RAND 

CARO  SFT  6  -  X  L  A  M 1  .  XL  AM? ,  X  I  n  c 

X L A M 1  =  LOWER  LIMIT  OF  THF  WAVELENGTH  FOR  THE  SPECTRAL 
X  L  A  m  R  =  IIDPPR  LIMIT  OF  THP  WAVFLFNGTH  FOR  THE  SPECTRAL 
XTMr  =  INCREMENT  SPECIFICATION  FOR  THF  WAVFLFNGTH 

CApn  SFT  7  -  LAMnni  (NOLAMS)  .  LAMRD2  ( NDL  AMS )  (FORMAT  3F10.fi) 

INTPGPATFO  RADIATION  LIMITS.  INPUT  ALL  LOWFR  THAN  ALL  UPPER 
VALUES.  THERE  WILL  RF  NOLAMS  VALUFS  FOR  FACH. 

capo  sft  r  -  temper  (napfas)  (format  bfio.6) 

TEMPERATURE  FOP  FACH  APEA.  INPUT  WOULD  CONSIST  OF  NIIME  CARDS 
WITH  NAPFAS  VALUES  PFP  CARD.  ^ 

CAPO  SET  9  -  (FORMAT  IS) 

BLANK  CAPO  WILL  FND  INPUT.  A  VALUE  TO  CONTINUE  READING, 
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Headings  are  provided  for  clarification  and  are  not  to  be  punched  data. 
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